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ABSTRACT 

We report the first ensemble results from the FUSE survey of molecular hydrogen 
in lines of sight with Ay ^ 1 mag. We have developed techniques for fitting computed 
profiles to the low- J lines of H2, and thus determining column densities for J = and 
J = 1, which contain >99% of the total H2. From these column densities and ancillary 
data we have derived the total H2 column densities, hydrogen molecular fractions, and 
kinetic temperatures for 23 lines of sight. This is the first significant sample of molecular 
hydrogen column densities of ~10^^ cm~'^, measured through UV absorption bands. 
We have also compiled a set of extinction data for these lines of sight, which sample 
a wide range of environments. We have searched for correlations of our H2-related 
quantities with previously published column densities of other molecules and extinction 
parameters. We find strong correlations between H2 and molecules such as CH, CN, and 
CO, in general agreement with predictions of chemical models. We also find the expected 
correlations between hydrogen molecular fraction and various density indicators such as 
kinetic temperature, CN abundance, the steepness of the far-UV extinction rise, and the 
width of the 2175 A bump. Despite the relatively large molecular fractions, we do not 
see the values greater than 0.8 expected in translucent clouds. With the exception of a 
few lines of sight, we see little evidence for the presence of individual translucent clouds 
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in our sample. We conclude that most of the lines of sight are actually composed of two 
or more diffuse clouds similar to those found toward targets like C Oph. We suggest a 
modification in terminology to distinguish between a "translucent line of sight" and a 
"translucent cloud." 

Subject headings: ISM: abundances — ISM: clouds — ISM: lines and clouds — ISM: 
molecules — ultraviolet: ISM 

1. Introduction 

Molecular hydrogen is the most abundant molecule in the interstellar medium, dominating 
the composition of the dense clouds that contain most of the mass. Even in diffuse clouds, H2 is 
present, containing from 10^^ to about half of the total hydrogen nuclei (e.g., Spitzer, Cochran, & 
Hirshfeld 1974; Shull & Beckwith 1982). Clearly, a full understanding of the physics and chemistry 
of the ISM requires a detailed knowledge of molecular hydrogen. 

As a homonuclear molecule, H2 has no dipole moment and hence no allowed rotational or 
vibrational transitions in the radio and infrared spectral regions. With the exception of forbidden 
quadrupole transitions in the infrared, which can be seen in emission in regions heated by radiation 
or shocks (Timmermann et al. 1996), or, in rare cases, in absorption when sufficiently high column 
densities are probed (Lacy et al. 1994), the only way to observe cold interstellar H2 is through its 
electronic transitions in the far ultraviolet. 

The Lyman (B^X) and Werner (C^X) bands lie in the spectral region between about 844 A 
and 1126 A. The moment of inertia for this low-mass molecule is small, resulting in widely separated 
rotational lines which are easily resolved. As a result, the far-UV spectrum of any reddened star is 
dominated by a wealth of H2 absorption bands. For a summary of the spectra of these bands, see 
Morton & Dinerstein (1976) or Barnstedt et al. (2000). 

Previous far-UV observations of H2 absorption have been conducted by various short-term 
missions, including sounding rockets (Carruthers 1970 - the first detection of H2 in space); the 
Hopkins Ultraviolet Telescope (e.g., Blair, Long, Sz Raymond 1996); ORFEUS (e.g., Barnstedt et 
al. 2000); and IMAPS (e.g., Jenkins &: Peimbert 1997). But by far the most extensive previous 
observations of far-UV H2 absorption were performed by the Copernicus satellite, which provided 
the first general quantitative studies of interstellar H2 as well as a wealth of information on its 
formation, its abundance, and its excitation in space (e.g., Spitzer et al. 1973; Spitzer, Cochran, & 
Hirshfeld 1974; Spitzer & Zweibel 1974; Jura 1975a,b; for summaries, see Spitzer k, Jenkins 1975 
and Shull & Beckwith 1982). The Copernicus work culminated in a survey of molecular and atomic 
hydrogen for 109 stars (Savage et al. 1977; Bohlin, Savage, &; Drake 1978). 

The Far Ultraviolet Spectroscopic Explorer {FUSE) instrument is well suited for observations 
of the H2 bands in absorption, and studies of H2 have been a longstanding goal of the FUSE mission 
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(Moos et al. 2000). FUSE covers the wavelength region from the Lyman Hmit at 912 A to about 
1187 A with spectral resolving power of about A/AA ~ 20,000. The throughput of FUSE has been 
specifically maximized in the middle of the Lyman band (around 1050 A), providing by far the 
most sensitive instrument yet available for observing cold H2 in the ISM. For a full description of 
the FUSE mission and instrumental performance, see Moos et al. (2000) and Sahnow et al. (2000), 
respectively. 

A subset of the FUSE Science team is conducting a program to study H2 in the densest clouds 
accessible, the so-called "translucent clouds" (van Dishoeck &; Black 1988). These clouds are defined 
as being optically thick (with total extinctions Ay > 1 mag) yet still sufficiently transparent as 
to allow optical and UV absorption-line measurements of interstellar atoms, ions, and molecules. 
The upper limit on extinction for a translucent cloud is usually considered to be about Ay ~ 5 
mag; beyond that limit there is usually not enough flux to allow high-resolution visible-wavelength 
spectroscopy, much less ultraviolet spectroscopy, which is always hindered by the rise in dust 
extinction at short wavelengths. Hence the goal of our program has been to observe and survey 
molecular hydrogen in translucent clouds, going as far as possible in the direction of maximum 
extinction, with the hope of penetrating clouds with Ay as high as 5 magnitudes. In doing so, we 
planned to extend the Copernicus-hased surveys of Savage et al. (1977) and Bohlin et al. (1978) 
to greater extinctions and denser clouds. 

More specifically, the goals of the FUSE translucent cloud H2 survey include: 

• Measuring total gas column densities, to help in determining interstellar gas-phase depletions 
and chemistry. 

• Studying the relationship between H2 and dust extinction, as an aid in assessing H2 formation 
models as well as extending correlations between extinction and H2 column density, sometimes used 
to estimate cloud masses. 

• Assessing the molecular fraction /h2 = 2A^(H2)/[2A/'(H2) -|- A^(H I)] and its relationship to 
dust extinction and other line-of-sight characteristics. 

• Measuring gas kinetic temperatures (Tqi) from the ratio of J = 1 to J = rotational states, 
and assessing the corrrelation of Tqi with other interstellar quantities. 

• Extending direct measurements of the CO/H2 correlation through UV absorption features 
of CO. This correlation is widely used to assess total masses of molecular clouds, but is currently 
based largely on CO abundances derived from mm- wave radio observations which do not necessarily 
sample the same material as H2 in absorption, and indirect H2 abundances estimated from dust 
extinction measures. 

• Assessing cloud physical conditions such as density and radiation field intensity by analyzing 
the excitation of high rotational states (J = 2 and greater). 

• Comparing H2 absorption measures from FUSE with infrared emission often attributed to 
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complex molecules (e.g., PAHs) as measured by the IRAS mission. 

• Measuring the abundance of HD in order to infer information on the D /H ratio. 

As will be seen later in this paper, all but the last three goals have now been met. We note 
also that UV absorption-line CO data do not yet exist for all of our program stars, so more progress 
on the CO/H2 correlation will come later. 

Two papers on FUSE observations of H2 in translucent clouds have already been published 
(Snow et al. 2000 [Paper I]; Rachford et al. 2001 [Paper II]). In addition to our translucent cloud 

H2 survey, another subset of the FUSE Science Team is conducting a general survey of molecular 
hydrogen in diffuse lines of sight, including Galactic stars with relatively little reddening, stars at 
high Galactic latitude, Magellanic Cloud stars (Tumlinson et al. 2002), and H2 in the lines of sight 
toward other extragalactic sources such as AGNs (Shull et al. 2000). 

The rest of the paper is organized as follows: in § 2 we describe the criteria for selecting target 
stars for the FUSE H2 survey, along with comments on the stars chosen; in § 3 we give details of 
the observations and data reduction and analysis; in § 4 we summarize the results; and in § 5 we 
discuss these results and their implications. 

2. Target Selection and Stellar Properties 

2.1. General Criteria for Selecting Target Stars 

Unlike the case of the diffuse H2 survey, which has been able to use spectra obtained by FUSE 
for other programs (e.g., the survey of stellar winds in OB stars; the survey of hot gas in the 
Galactic halo; and the survey of interstellar gas at high latitudes and along lines of sight toward 
the Magellanic Clouds and other galaxies), our sole source of data comes from stars we observe 
specifically for our program on translucent cloud lines of sight. No one observes heavily reddened 
stars unless absolutely necessary, because the dust extinction acts to increase the observing time 
needed, particularly at short wavelengths. Thus, wc had to conserve observing time and be as 
careful as possible in choosing our target stars while seeking to fulfill our basic criteria. 

The result was a list of 45 stars; 31 assigned to FUSE program P116, 4 assigned to QlOl, 
9 assig ned to P216i\ and one target (HD 24534) added from the personal program of one of us 
(T.P.S.; P193). Of these 45 stars, 25 had been observed through June 2001, and 23 are included in 
this paper. 



As a historical aside, the P216 program came about once FUSE had been in operation for several months, 
and it was established that the observing efEciency was better than expected. PI Team members were afforded an 
opportunity to request more guaranteed time for targets that did not conflict with other programs, and we selected 
the additional 9 targets, two of which are included in the present work. 
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The selection of target stars for the FUSE translucent cloud H2 survey program was based on 
several criteria: 

• Maximizing the total extinction probed; 

• Exploring lines of sight known to have a wide range of extinction characteristics such as Ry 
(the ratio of total to selective extinction) and the extinction parameters defined by Fitzpatrick & 
Massa (1986, 1988, 1990) on the basis of lUE data; 

• Relative simplicity of line-of-sight velocity structure based on high resolution observations 
of K I, Na I, and CH; 

• Availability of (or feasibility of obtaining) ancillary data on optical interstellar lines of atomic 
and molecular species not observed in the FUSE passband (including ultra-high resolution spectra 
for velocity structure analysis), and the availability (or feasibility of obtaining) infrared spectra of 
both gas-phase and solid-state absorption features. 

The first of these criteria was the most critical, due to the trade-off between maximizing 
dust extinction and finding stars with sufficient far-UV flux to be observable with FUSE within 
reasonable exposure times. Thus, we explored known UV fluxes for candidate stars, in most 
cases having to extrapolate to shorter wavelengths from the lUE cutoff at about 1170 A. The 
extrapolations were based on a combination of the known flux in the lUE band, the spectral type 
and instrinsic flux distribution of the star, and the shape of the UV extinction curve. 

In our initial selection of targets we drew heavily upon the lUE-based ultraviolet extinction 
curve survey of Fitzpatrick &; Massa (1986, 1988, 1990) and several papers on optical interstellar 
absorption lines, both atomic and molecular (e.g., van Dishoeck & Black 1989; Gredel, van Dishoeck, 
& Black 1993). Wc eliminated stars known to have complex line-of-sight Dopplcr structure due to 
multiple clouds. However, we found that this distinction was almost futile, as nearly every line of 
sight turned out to have complex structure when examined at sufficiently high spectral resolution. 
For example, the star we chose to analyze for our first paper from this program, HD 73882 (Paper 
I), turned out to have no fewer than 21 identifiable Doppler components in the ultra-high resolution 
spectra of optical interstellar lines due to K I and Na I. 

One of us (D.E.W.) has led a program to obtain very high-resolution optical spectra of our 
FUSE candidates, and has succeeded in obtaining spectra at resolving powers of A/ A A > 150,000 
in most cases. The results will appear in a separate paper (D. E. Welty et al. 2002, in preparation) 
while being made available to us as we analyze the FUSE H2 spectra. But, these high-resolution 
spectra are not very important for the current survey, since we include here only the two lowest 
rotational states of H2, whose absorption lines are damped and which therefore can be analyzed 
accurately without concern for the detailed velocity structure. We do, however, reference the optical 
spectra in a few cases. 

Table 2 provides ancillary data, where available, on carbon-based molecular species for the 
stars listed in Table 1. Additional data regarding the CH component structure will be discussed 
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in § 5.1. We plan to supplement Table 2 by performing our own ground-based observations of 

various atomic and molecular species. We expect especially to emphasize infrared spectroscopy in 
this effort, in order to obtain data on both gas-phase and solid-state absorption features such as 
those due to CO, water ice, and the 3.4-//m hydrocarbon feature. 

2.2. Extinction properties 

The fundamental quantities used to categorize dust extinction along a line of sight are E{B — V) 
and Rv = Av/E{B — V), where Ay is the total visual extinction in magnitudes. Since E{B — V) 
is also expressed in magnitudes, Ry is unitless. The quantity E{B — V), or color excess, represents 
th(^ amount of reddening along a line of sight and gives a measure of the amount of dust and gas. 
The quantity Ry, or total-to-sclcctivc extinction ratio, gives a measure of the size of dust grains, 
and is a convenient single parameter for describing the overall shape of extinction curves (Cardelli, 
Clayton, & Mathis 1989). The detailed shape of the UV extinction curve also carries considerable 
information. We discuss these extinction properties in the following sections. 

2.2.1. E{B-V) 

The determination of the E{B — V) color excess is relatively straightforward for early-type 
stars. The intrinsic {B — V)q colors for O and B stars vary slowly as a function of spectral type, 
minimizing errors due to inaccurate typing and variations in Ry. Table 3 gives our adopted values. 
In most cases, these values come from either an analysis of the upper main sequences of OB 
associations, or simply from a comparison of the observed colors to tabulated intrinsic colors for a 
given spectral type. The dift'crcnces between the various references for a given star are typically only 
a few hundredths of a magnitude at worst. This is comparable to the differences between different 
lists of intrinsic colors, and not much greater than the photometric uncertainties themselves. In the 
case of the Be star HD 110432, the value quoted in Table 3 includes a correction to the observed 
E{B — V) to allow for the emission line behavior. 

2.2.2. Ry 

The situation with regards to Ry is not as simple. Whatever the method, measuring Ry 
requires more difficult observations with much greater calibration problems. We will exploit three 
methods for estimating Ry in the present work. 

If a star does not show excess infrared emission (e.g. Be stars or other stars with circumstellar 
material), an extrapolation of infrared color excesses to infinite wavelength yields an estimate of 
Ry. In the method developed by Martin & Whittet (1990), the near-infrared extinction curve takes 
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the form, 

Ex-v 



Eb-v 

The infrared extinction in the range 0.9-4.8 /xm has been expressed as 



= eA-'* - Rv- (1) 



Eb-v 



= eA-", (2) 



where e is the power-law amphtude, similar to the C4 parameter of the UV extinction curves 
described in § 2.2.3, and A is given in fim.. Martin & Whittet found that the power-law index, a, 
rarely deviates from the average value of 1.8, even for material with Ry far from the interstellar 
average of 3.1. Thus, to within the normalization provided by e, all extinction curves in the range 
A = 0.9-4.8 //m can be described by this universal curve. With the known quantities Ex^v/Eb-v 
and A^'^ as the ordinate and abscissa, respectively, wc can fit a linear relation with slope e and 
a ?/- intercept of —Ry- This allows us to apply Eq. 1 to cases where we only have 3 infrared 
measurements (or even 2 measurements, with much greater uncertainty). We have verified that in 
the cases where we have 4 or 5 measurements, the full form of Eq. 1 gives results consistent with 
the linearized form; i.e. a 1.8. Although Eq. 1 and 2 appear to hold for photometric bands from 
/ through M, we have limited our fits to the five available bands from J through M. As Table 3 
indicates, we generally had 3 or 4 bands to work with. 

One complication of this method is that the colors of the stars must be corrected for reddening 
by comparison with a standard relationship as a function of spectral type. In absolute terms, the 
infrared colors (J — V)q, (H — V)q, etc. change much more rapidly with changing spectral type 
than (B — V)q. The choice of the appropriate set of intrinsic colors, as well as the spectral type, 
leads to significant differences in the derived value of Ry- Many authors have used older intrinsic 
color-spectral type relationships such as that of Johnson (1966). However, one must interpolate to 
generate the appropriate H — V colors, and the tables are incomplete for O stars. More recently, 
Wcgncr (1994) compiled a more complete set of intrinsic colors for O and B stars, and other similar 
calibrations exist. 

We found infrared photometry in at least 2 filters for 19 of our 23 stars and performed fits 
to Eq. 1 based on intrinsic colors from both Johnson (1966) and Wegner (1994). In many cases, 
photometric uncertainties for a specific star were not given. Thus, in our error analysis, we assumed 
a = 0.03 for JHK, a = 0.05 for LM, and a = OME{B - V) for E{B - V). The final statistical 
uncertainties on Ry were then close to 0.2 in all cases. 

In two cases, HD 110432 and HD 168076, we could not obtain reasonable fits to Eq. 1. HD 
110432 is a Be star, so that we would not expect this method to work due to the distorted infrared 
colors. HD 168076 also shows emission lines. However, several other O stars in our sample also 
show emission lines, yet the infrared photometry closely follows the expected power law. We have 
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noticed for this star that the infrared photometry from The et al. (1990) differs significantly from 
that given by Aiello et al. (1988), yet neither set of photometry obeys Eq. 1. 

For the remaining stars with photometry in 4 or 5 filters and good fits, the differences between 
the two sets of intrinsic colors are large and systematic. The differences averaged 0.20, and in all 
cases the values of Ry were smaller for the Wegner (1994) intrinsic colors. 

For the 7 cases where previous authors have used this or a similar method (BD +31° 643, HD 
62542, HD 73882, HD 206267, HD 207198, HD 210121, and HD 210839), our values agree nearly 
exactly when we have used the same set of intrinsic colors. Considering that the differences in Ry 
for the two different sets of intrinsic colors are smallest for the B stars, and that the Wegner (1994) 
calibration is much more complete for the O stars, we have accepted this calibration as the more 
accurate overall and give those values in Table 3. 

The wavelength of maximum linear polarization generally correlates well with Ry, following 
the relationship 

Ry = (5.6 ± 0.3)A^ax, (3) 

where Amax is given in ^um (Whittet & van Breda 1978). Combined with the uncertainties in 
determining Amax (a few tenths of a percent to several percent), this method produces values of Ry 
accurate to about 5-10%, similar to the uncertainties from the first method. 

It is worth emphasizing that the relationship between Ry and Amax was calibrated in part 
from infrared photometry, so that this second method is not completely independent of the first. 
A further complication is that Eq. 3 may not always be satisfied. Recently, Whittet et al. (2001) 
found deviations from Eq. 3 within the Taurus dark clouds for a number of lines of sight, in the 
sense of larger than normal Amax for normal values of Ry derived from IR photometry. 

Finally, one can estimate Ry by comparing the ultraviolet extinction curves to standard curves. 
Again, this method requires calibration from other methods of determining Ry. A comparison of 
the shapes of extinction curves can be somewhat subjective. However, in an important set of papers, 
Fitzpatrick & Massa (1986, 1988, 1990) devised a six-parameter description of the ultraviolet curves 
that could reproduce the entire curve to within the observational errors. We discuss extinction 
curves in more detail in the following section, but the most useful correlation is between Ry and 
the linear far-UV rise, the "02" parameter. The relationship derived by Fitzpatrick (1999) is 

C2 = -0.824 + 4.717ii:-^ (4) 

The typical scatter about this relationship is a few hundredths in Ry^, or ~0.2 in Ry for small 
values, and ~0.5 for large values. Larger deviations occur for a few lines of sight. For example, the 
extinction curve for HD 62542 has an extremely steep far-UV rise, out of proportion with the rest 
of the curve (Cardelli &; Savage 1988). 
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Unfortunately, in some cases, there are significant disagreements in Ry derived from the differ- 
ent methods. However, in most cases the agreement between the methods is reasonable. The most 
direct method for the determination of Ry is the method involving IR photometry, and we have 
such a determination for a majority of our targets. We have used this method as our primary source 
for Ry, using values from polarimetry and the extinction curves in that order to supplement the IR 
photometry. We then simply derive Ay from the product of Ry and E{B — V) (Table 3). We have 
not included specific uncertainties on these values, but as wc noted above, the typical statistical 
uncertainty on the photometric values of Ry is 0.2. The resultant statistical uncertainties on Ay 
are then typically 0.1-0.3. Systematic effects may be important as noted above in the discussion 
of the methods of determining Ry . 



2.2.3. Extinction curves 

As mentioned previously, Fitzpatrick & Massa (1986, 1988, 1990) devised a six-parameter 
scheme to describe the UV extinction curve in the lUE wavelength range, ~ 1170-3200 A. The six 
parameters include three that describe the central wavelength (xq), width (7), and height (03) of 
the 2175 A bump, two that represent a linear background term (ci and C2), and one that describes 
the far-UV curvature (04). With x = in units of /xm~^, the function describing the extinction 
curves, k, is given by 



k{x -V) = ci + C2X + C3D{x; 7, xq) + C4F{x), (5) 

where 



a;2 



Dix;j,xo)= (6) 



and 



, I 0.5392(x - 5.9)2 + 0.0564(x - 5.9)=^ x > 5.9 
10 X < 5.9 

The function D{x;^, xq) is the so-called Drude function (the absorption profile of a damped 
freely oscillating electron) and represents the 2175 A bump. The far-UV curvature is described 
by the function F{x). As previously noted, C2 correlates well with Ry, and grain size. Table 4 
gives previously published extinction curve parameters for most of our targets. In § 4.5 we will 
investigate various correlations involving these parameters. 
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2.3. Special line-of-sight characteristics 

Many of the lines of sight are of special interest due to their location and/or environment. 
We give a brief overview of each line of sight in the following sections, including the mention of 
particularly interesting values from Tables 1-4. 

2.3.1. BD +31° 643 

The binary system BD +31° 643 lies within the cluster IC 348 and interstellar material as- 
sociated with the cluster. The line of sight is only 8' from the well studied line of sight toward o 
Per. However, the latter star lies in the background of IC 348 and its associated material (Sancisi 
et al. 1974). Snow ct al. (1994) found a "composite" extinction curve, with contributions from 
material having a broad 2175 A bump and steep far-UV rise, and from material having a narrow 
bump and shallow far-UV rise. Also, the authors noted very different molecular abundances in the 
two lines of sight: enhanced CH and CH+ column densities toward BD +31° 643, and a reduced 
CN abundance. They concluded that star formation in this cluster has altered the physical and 
chemical state of the gas. Snow (1993) has argued that the enhanced CH+ column density is larger 
than would be expected from shock formation, and proposed a radiative formation mechanism due 
to a very strong UV radiation field. 

BD +31° 643 is surrounded by a circumstellar disk, only the second such disk ever directly 
imaged (Kalas & Jcwitt 1997). Andersson & Wannier (2000) suggest that more then one magnitude 
(approximately half) of the total visual extinction may arise in this disk. This might lead to a 
composite extinction curve, and a single Ry might not be an appropriate description of the line of 
sight. 

2.3.2. HD 24534 

This Be star is more commonly known as X Persei and has an X-ray-bright pulsar companion. 
However, the line of sight toward the pair intersects a molecular cloud, and the UV flux is among 
the largest known for a target with Ay > 1 mag. In fact, it was (barely) bright enough to be 
observed with Copernicus (Mason et al. 1976). They found log A''(H2) = 21.04, the largest such 
value obtained with Copernicus. The molecular fraction was also estimated to be very large (Snow 
et al. 1998), albeit the H2 column density was rather uncertain. 

2.3.3. HD 27778 

This binary star lies within or just behind the Taurus- Auriga molecular cloud complex (Kenyon, 
Dobrzycka, & Hartmann 1994), sampling relatively low-extinction material in this complex (Whittet 
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et al. 2001). Although fit parameters have not been pubUshed for the extinction curve, the curve 
depicted in Joseph et al. (1986) indicates greater-than-normal far-UV extinction, and a shallower- 
than-normal 2175 A bump, characteristic of a "dense cloud" environment. 

We note that while Joseph et al. (1986) reported an H I column density log N = 21.59±0.11 
from measurement of the Lyct line, at spectral type B3 V the Lyman-series interstellar lines will be 
strongly contaminated by stellar lines (Savage &: Panek 1974; Shull &: van Steenberg 1985; Diplas 
& Savage 1994). 

2.3.4. HD 62542 

HD 62542 lies in or just beyond the Gum nebula, and the line of sight passes through an isolated 
ridge of molecular material (Cardelli et al. 1990; O'Donnell, Cardelli, &; Churchwell 1992). Despite 
a normal value of Ry, Cardelli &; Savage (1988) found that the UV extinction curve is peculiar 
in two important respects. The far-UV portion of the curve is unusually steep, and the 2175 A 
bump is weak and strongly shifted to 2110 A. Combined with the following additional observations, 
Cardelli et al. (1990) concluded that the line of sight contains very high density material. First, 
optical maps of the Gum nebula complex indicate small, dark globules scattered through the ridge of 
material through which the light from HD 62542 passes. Second, the CN excitation temperature is 
considerably larger than the temperature of the cosmic background radiation, indicating significant 
collisional excitation. Finally, the column densities for CN and CH are exceptionally large given 
the relatively small visual extinction. The conclusion reached by Cardelli et al. (1990) was a 
remarkably high density of n ~ 10^ cm~^, and a cloud with thickness d ~ 0.2 pc and mass m ~ 1 

Other density diagnostics suggest less extreme conditions, i.e., C2 excitation implies n ~ 500- 
1000 cm~^ (Gredel et al. 1991, 1993), and CN observations also imply n ~ 500 cm~^ (Black & van 
Dishoeck 1991). With this value, the cloud size is larger, i.e., 0.5-1.0 pc (Snow et al. 2002). In any 
case, the presence of such high densities in a cloud with relatively small extinction is unusual. 

Despite the large column densities of CN and CH, the CH+ column density is remarkably 
small. If shocks are the main formation mechanism for CH+, this result would be consistent with 
high-density material. However, as discussed in § 4.1 the formation mechanism of CH+ is still in 
question. 

2.3.5. HD 73882 

Molecular hydrogen has been studied in some detail for this line of sight in Paper I. We take 
advantage of additional data obtained subsequent to this paper, but this information does not alter 
the original conclusions. The extinction is among the largest in the present study. The extinction 
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curve shows a shallow 2175 A bump with steeper than normal far-UV extinction, characteristics 
common to "dense cloud" curves. However, the overall conclusion of Paper I was that the line of 
sight appeared to be similar to diffuse clouds, and may contain a collection of several clouds. 

2.3.6. HD 96675, HD 102065, HD 108927 

These three targets all lie within or just beyond the Chameleon dark clouds. Although the 
color excesses are among the smallest of the present sample, the column density of the CH radical 
toward HD 96675 and HD 102065 suggests a considerable molecular content (Gry et al. 1998). 
High- velocity atomic gas is also seen toward these two stars (Gry et al. 1998). Gry et al. (2002) 
have studied H2 toward these lines of sight and concluded that a component of the high- J excitation 
comes from coUisional excitation in regions of warm gas. 

Despite their relative proximity, the three lines of sight sample material with strongly differing 
dust properties, as indicated by IRAS photometry, the UV extinction curves, and Ry (Boulanger 
et al. 1994). Unfortunately, the methods for the determination of Ry given in Table 3 show 
exceptional disagreement. 

2.3.7 HD 110432 

This line of sight has been studied in detail in Paper II. The main conclusion was that the 
line of sight was similar to C Oph, based on the column densities in the J = 0-7 levels, as well as 
chemical modeling. The CN column density is quite small relative to our present sample, but is 
comparable to that seen toward ( Oph. 

2.3.8. HD 154368 

This star is one of the few with sufficient flux to be observed by FUSE through 2.5 magnitudes 

of visual extinction. The component structure appears to be dominated two components, and many 
species have been well studied with HST (Snow et al. 1996). The extinction curve is more or less 
normal, but tends slightly towards the "dense cloud" curves (Snow et al. 1996). They concluded 
that the line of sight contains extended regions of moderately dense gas, as opposed to one or more 
dense cloud cores. 

2.3.9. HD 167971 

This target is actually an eclipsing binary pair, with a third supergiant component which is 
the most luminous member of the group (Leitherer et al. 1987). The three stars may represent 
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a Trapezium-type system or possibly a chance alignment. All three stars are likely members of 
the Serpens 0B2 association, and at least the supergiant star lies within the cluster NGC 6604. 
Serpens 0B2 is physically related to the H II region Sharpless 54. This target has the largest color 
excess and visual extinction of the present sample. 

2.3.10. HD 168076 

This star lies within the cluster NGC 6611, which in turn is associated with the dusty H II 
region M16. At type 04, this star is the earliest target included in our survey. The cluster stars 
in NGC 6611 show significant intracluster and/or foreground differential reddening (Orsatti, Vega, 
&; Marraco 2000). These authors found evidence for a slightly larger value of the wavelength of 
maximum polarization, \max^ 

than in the general ISM. However, the line of sight toward HD 168076 
did not conclusively show this effect. 

2.3.11. HD 170740 

Although a bright target in both the optical and ultraviolet, there have been relatively few 
studies of this line of sight. The column densities reported in Table 2 are unremarkable and the 
UV extinction curve was judged normal by Witt, Bohlin, & Stecher (1984). 

2.3.12. HD 185418 

This line of sight is not well studied. Fitzpatrick & Massa (1986, 1988, 1990) included the 
target in their extinction curve survey. The extinction curve shows a stronger than normal 2175 A 
bump, along with less than normal far-UV extinction. 

2.3.13. HD 192639 

HD 192639 lies within the cluster NGC 6913. Nichols-Bohhn & Fesen (1993) studied the 
complex interstellar environment centered on the Wolf-Rayet star HD 192163, lying just one degree 
away from HD 192639 and at nearly the same distance. The authors found evidence for a supernova 
shell surrounding HD 192163, and perhaps HD 192639. Most of the lines of sight toward hot 
stars in the vicinity showed high-velocity components (primarily at negative velocities) including 
HD 192639, indicative of the multiple superbubble structure surrounding the Cyg OBI and 0B3 
associations. 

P. Sonnentrucker et al. (2002; in preparation) have undertaken a detailed study of abundances 
along this line of sight, combining the FUSE data with high-resolution HST spectroscopy. Based 
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on the inferred low particle density, complex component structure, and chemical depletions, they 
concluded that the line of sight is a collection of diffuse clouds. 

2.3.14. HD 197512 

This is another line of sight from the Fitzpatrick and Massa (1986, 1988, 1990) sample that 
otherwise has not been well studied. The extinction curve shows greater extinction than the 
standard curve at all wavelengths, suggesting a possible miscalibration with E{B — V) or spectral 
type of the comparison target. 

2.3.15. HD 199579 

This is the second target in the present survey for which A'^(H2) was also derived from Coperni- 
cus observations. The star contributes to the excitation of the North America Nebula (NGC 7000), 
but is not the primary excitation source (Neckel, Harris, & Eiroa 1980). The extinction curve 
shows larger than normal far-UV extinction, but the 2175 A bump is normal. The line-of-sight CN 
column density is unusually small. 

2.3.16. HD 203938 

The extinction curve for this line of sight is similar to that for HD 197512 in the sense that 
the extinction is greater than normal at all wavelengths (Fitzpatrick &: Massa 1990). However, the 
effect is not as great for this target. Otherwise, the line of sight is not well studied. 

2.3.17. HD 206267 

HD 206267 is a Trapezium-like quadruple star system (Abt 1986) within the cluster Trumpler 
37. The cluster is associated with the H II region IC 1396, and the hottest star of the group, HD 
206267A, is the main exciting source for the region. The FUSE observation contains contributions 
from components A and B. The cluster and H II region have been well studied. Morbidelli et 
al. (1997) found relatively uniform extinction across the cluster both in terms of total extinction 
and the total-to-selective extinction ratio. Clayton & Fitzpatrick (1987) derived UV extinction 
curves and found that HD 206267AB, like most early-type members of the cluster, shows a normal 
2175 A bump and stronger than normal far-UV extinction. The uniformity of extinction and 
extinction curves for the cluster stars led to the interpretation that most of the extinction toward 
HD 206267AB is foreground with a small contribution from IC 1396 (Clayton k, Fitzpatrick 1987; 
Morbidelh et al. 1997). 
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More recently, however, variations in the column densities of individual velocity components 

over spatial scales of ~10, 000-40, 000 AU have been observed toward the HD 206267 system (Lau- 
roesch &: Meyer 1999; Pan, Fcdcrman, Welty 2001). The most striking variations appeared for CN, 
the most density sensitive of the species studied. 

2.3.18. HD 207198 

The extinction curve derived by Jenniskens &: Greenberg (1993) indicates much stronger than 
normal far-UV extinction, and the other indicators of Ry also indicate that this is an unusual line 
of sight. In this regard, it is similar to the line of sight toward HD 210121, but less extreme (see 
below) . 

2.3.19. HD 207538 

The polarization data suggest that this line of sight may be similar to HD 207198 and HD 
210121, with a small value of Ry- However, the line of sight is otherwise poorly studied, and we 
are unaware of a published extinction curve that might verify the unusual extinction. 

2.3.20. HD 210121 

HD 210121 lies within or behind the remarkable high-latitude molecular cloud DBB 80 (Desert, 
Bazell, & Boulanger 1988; de Vries & van Dishoeck 1988), and this is the only high-latitude line of 
sight in the present survey. The UV extinction curve exhibits one of the steepest far-UV rises ever 
observed (Welty & Fowler 1992), consistent with the exceptionally small total-to-selective extinction 
ratio, Rv = 2.1. In addition, the 2175 A bump is very weak, such that all available extinction 
parameters suggest dense material. This line of sight is especially important as it is the only one 
known with such extreme extinction that is bright enough for FC/S'i? observations. While HD 62542 
has a similar far-UV rise, it does not have the corresponding extreme Ry- Welty & Fowler (1992) 
also found a small ratio between the 100 /xm flux and the total hydrogen column density, suggesting 
a smaller than average incident radiation field consistent with the cloud's location 150 pc from the 
Galactic plane. 

Grain models suggest an excess of small grains and a deficiency of large grains relative to 
the average interstellar medium (Larson et al. 2000). Li & Grccnburg (1998) have modeled the 
extinction, polarization, and emission in the molecular cloud by including a contribution from grains 
which have undergone erosion and thus have thinner mantles. 
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2.3.21. HD 210839 

Also known as A Cephei, this is the third and final target in our list for which iV(H2) was also 
derived from Copernicus observations. As the optically brightest target in the present survey, many 
ground-based studies have been performed. However, the column densities of the various species 
(Table 2) are generally normal. Bless k. Savage (1972) published a partial UV extinction curve out 
to 1800 A which indicates a nearly normal 2175 A bump. 

Jenkins &; Tripp (2001) found two main components in high resolution observations of C I 
lines (with possible structure within each component). The two components were separated by 
about 23 km s~^. The weaker red component, presumably associated with the expanding Cepheus 
bubble, contains high-pressure gas {p/k = 10^'^ cm~^ K), as derived from C I and O I fine-structure 
excitation. This component would be a good candidate for the presence of H2. While the 23 km 

separation could be resolved in weak H2 lines in our FUSE spectra, saturation broadening in 
lines with J < 4 (and the strongest lines with J = 5) overwhelms any potential weak red-shifted 
component. In the weaker lines with J = 5-7 we do not see resolved component structure. Based 
on the strong CH (and CH^) absorption seen by Crane, Lambert, & Sheffer (1995) at the same 
velocity as the blue components of C I and O 1 seen by Jenkins & Tripp (2001), we assume that all 
the observed H2 is also associated with this component. 

3. Observations and data analysis 

Table 5 gives information on our FUSE observations. This survey is based on a total of 288 
ksec of integration time on 23 targets, taken over the course of 20 months. Our targets are generally 
relatively faint, and only HD 24534, HD 110432, HD 170740, HD 199579, and HD 210839 were 
bright enough to require spectral-image mode (Moos et al. 2000). 

BD-|-31° 643 and HD 73882 were observed during the "early release observation" phase of the 
mission, and were subsequently reobservcd due to incomplete data coverage. The first half of the 
initial observation of HD 207538 was lost due to a software problem onboard the spacecraft, and 
this target was also reobserved. In several cases, the steepness of the extinction curve and/or the 
faintness of the target prevents us from having usable data in the shortest wavelength segments 
(SiC IB and 2A). However, the wavelengths below 980 A are not used in the present H2 analysis. 

As we are fitting very broad profiles and our earliest datasets have poor S/N, we have not 
made an attempt to re-process all of the data sets with the most recent version of the CALFUSE 
pipeline. We have, however, applied the most appropriate version of the wavelength calibration to 
each dataset. 

In all cases, each observation is broken into two or more individual integrations. Before com- 
bining the individual spectra, we perform a cross-correlation analysis on a cluster of narrow lines 
near the center of each data segment to co-align the spectra. Since the spectra are highly over- 
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sampled (one resolution element corresponds to about 9 pixels), we only applied integral pixel 
shifts, so that the data are never interpolated or resampled. The pixel-by-pixel uncertainties are 
propagated through the co-additions, and each of the 8 detector segments is processed separately. 
Significant differences in flux calibration, line spread function, and wavelength solutions conspire 
to make useful co-additions of data from different segments difficult. Despite these differences, we 
do find generally good agreement between segments; however, the differences are larger than the 
formal uncertainties within each segment, as discussed in the following section. 

We have reported per-pixel values of S/N in Table 5. With the typical 9 pixel resolution 
element, the S/N for one resolution element could be 3 times as large. The true S/N for a resolution 
element for our highest quality observations may be limited by fixed-pattern noise and the lack of 
a flat-fielding correction. However, thermal motions present in the mirrors and gratings help to 
dither the spectral image across the detectors over the course of multiple integrations and help 
smooth out the fixed-pattern noise. Thus, in most cases our estimated S/N for a resolution element 
is close to the optimum value. More information concerning the on-orbit performance of FUSE is 
given by Sahnow et al. (2000). 

3.1. Data Analysis 

3.1.1. Overview of the problem 

A description of our analysis procedures is given in Paper II. However, in the present paper 
we expand upon several issues involved with fitting these complex spectra. This discussion is 
based on a long period of experimentation with various analysis techniques on spectra with a wide 
range of S/N and column densities. These techniques were also used by Tumlinson et al. (2002), 
with appropriate modifications to allow for the presence of both Galactic and Magellanic Cloud 
components in the overall H2 spectrum. 

Figure 1 shows sample spectra for a high S/N target (HD 210839) and a low S/N target (HD 
154368). The equivalent widths of the undamped J > 3 lines can generally be measured individually 
without regard to the specifics of the unresolved component structure along the line of sight, and 
then a curve of growth analysis can be performed. We are performing such an analysis, and these 
results will be presented at a later date. However, to determine the total H2 column density we 
must perform profile fits on the highly damped J = and 1 profiles. In addition, since the R(2) 
lines are blended with the main J = and 1 profiles, we must include J = 2 in the profile fits. 
Fortunately, the P(2) lines are sufficiently isolated to constrain the J = 2 column densities for our 
present purposes. 

The J = and 1 lines themselves are too heavily damped to be sensitive to the detailed 
component structure (or 6- value if a single absorbing cloud is assumed). The J = 2 lines are 
somewhat sensitive to this choice, and in turn the blending between the P(l) and R(2) lines can 
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afFect the derived J = 1 column density. However, we have found that a change in 6- value of a 
factor of 2 typically changes log iV(l) by just 0.01 dex, even with a much greater change in N{2). 

The steep extinction curves in these heavily reddened lines of sight and the overlap between 
Lyman and Werner bands of H2 and the Lyman H I lines prevent us from reliably using H2 bands 
near and shortward of Ly/3. Thus, we have considered Lyman bands of H2 from (0,0) through (4,0). 

The continuum radiation from the background hot stars is punctuated by photosphcric lines. 
If these lines lie in the far wings of the damped profiles, they can be easily divided out of the 
spectra. If they lie near the zero-intensity cores of the profiles, these lines do not affect the overall 
H2 spectrum. However, contamination of portions of the H2 profiles with normalized intensities of 
~0.5 is both difficult to remove and can cause significant disruption of the H2 spectrum. In Paper 
II we noted that fits of the Lyman (0,0) band in HD 110432 exhibited such contamination, and 
that the (3,0) band fits were also possibly affected. Further exploration of this issue in additional 
spectra indicate that contamination of these bands is common. While detailed modeling of the 
stellar spectra is beyond the scope of this work, we note that the (1,0), (2,0), and (4,0) bands 
appear to be the cleanest of the long wavelength bands in hot star spectra observed through less 
H2. In the interest of producing a uniform measurement of H2 column densities across our sample, 
we are limiting our J = and 1 analysis to these bands. These bands appear a total of nine times 
on five different detector segments. 

A final problem is the wide range in data quality. As seen in Table 5, the per-pixel S/N varies 
from unity to nearly 30. Our data analysis techniques have to work well across this large range. 

3.1.2. Fitting techniques 

Our goal is to match a model spectrum of the low-J lines to the data by varying the J = 
0, 1, and 2 column densities, a quadratic continuum, and a zero-point wavelength shift. We have 
applied two distinct fitting methods to our data to minimize the squared difference between the 
model and the data, non-linear least squares (i.e., the Levenburg-Marquardt "CURFIT" algorithm 
from Bevington & Robinson 1992), and the "downhill simplex method" (the "AMOEBA" algorithm 
from Press et al. 2000). The non- linear least squares method has the advantage of producing formal 
uncertainties on each fit parameter from the covariance matrix when using the appropriate weighting 
scheme. In addition, it is much less computationally intensive. However, this method also requires 
the evaluation of partial derivatives with respect to each parameter, but our modeled profiles do 
not have analytical derivatives. The downhill simplex method works to minimize a quantity, in 
this case the difference between the model and the data, weighted in some manner. The method is 
more computationally intensive but only requires function evaluations and not derivatives. It can 
also be much more robust, particularly when dealing with a large number of fit parameters. 

We have experimented with several weighting schemes for the data, including "instrumental" 
(1/error^), "signal-to-noise" (data/error), and "uniform". Instrumental weighting has the advan- 
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tage of producing genuine values and appropriate statistical error bars on each fit parameter, 
and is the usual choice for astronomical data. However, we have many data sets with very poor data 
quality (per-pixel S/N of a few or less). In these cases, the generally Poissonian errors can not be 
reasonably approximated by a normal distribution. This fact tends to skew the weighted fits toward 
the pixels with smaller values instead of bisecting the data points. Thus, in the poorest quality 
observations in the present sample, the instrumentally weighted fits do not provide a good match 
to the data. If we re-bin the data to increase the S/N, this becomes less of a problem. However, 
such re-binning has its own set of problems. Signal-to- noise weighting is somewhat less susceptible 
to this problem, although this weighting method is rarely used. A uniformly weighted fit removes 
this effect completely, but gives limited information on the uncertainties of the fit parameters. In 
the uniformly weighted fit we are in effect assuming that the perfect model should bisect the data 
points in a given range for any S/N. 

We have tested the various combinations of fitting techniques and weighting schemes. As an 
example, in Table 6 we show the results of a suite of fits of the (1,0), (2,0), and (4,0) bands from 
the LiF lA and LiF 2A segments for HD 199579. This is a high-quality spectrum with S/N of 20 

per-pixel at the maxima between the H2 bandheads, and we would expect to see good agreement 
among the various methods. The values in Table 6 show that this expectation is verified. Thus, we 
can choose the fitting technique and weighting scheme most appropriate for our poor-quality data 
without sacrificing the accuracy of the fits to the good-quality data. 

As we observed in Paper II, the largest variation in the fit parameters occurs when comparing 
one band to another. In addition, the formal uncertainties we derive for the high-quality data sets 
through both the CURFIT method are less than 0.01 dex. We also performed Monte Carlo simu- 
lations whereby we either added noise to synthetic profiles matching the data, or added additional 
noise to the data itself, and again the changes in the column densities were much smaller than the 
band-to-band and segment-to-segnient differences. From these findings, we confirm our conclusion 
in Paper II that the largest source of errors in the low-J column densities are from effects other 
than Poisson noise. These effects include low-level contamination from stellar lines, fixed-pattern 
noise, wavelength-dependent errors in the flux calibration, and other factors. 

3.1.3. Summary of the fitting technique 

Based on the previous discussion we now describe our revised fitting technique in full. Before 
fitting the spectra we identify all visible atomic lines and all visible H2 lines with J > 3, model 
them with Gaussian profiles, divide the lines out of the spectrum, and exclude from further fitting 
the cores of the removed lines if they dropped below half of the local continuum level. We also 
remove obvious stellar lines with a similar procedure. This leaves only seven fit parameters: a 
quadratic polynomial continuum, logarithmic column densities for J = 0, 1 and 2, and a zero-point 
wavelength shift. In Paper II, we included the high-J and atomic lines in our fits instead of an 
outright removal of the lines. Our choice in this matter has little effect on the quality of the fits, but 
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the reduction in fit parameters if we remove the lines reduces the computational time, and tends to 
make CURFIT as robust as AMOEBA. Also, the number of individual lines that must be modeled 
is reduced, again reducing the computational time. During the line-removal phase, we also select 
the appropriate wavelength ranges for the fitting of each band. We make these ranges as uniform 
as possible from target-to-target and for fits of the same band in different detector segments of the 
same observation. 

We assign a single 6-value to represent the overall component structure, which only affects the 
modeling of the J = 2 lines. As stated previously, this value need not be very accurate, and we 
have used preliminary values from our curve-of-growth analysis of the high- J lines, combined with 
the high-resolution ground-based data described in § 1. 

While there must be unresolved velocity structure in nearly all cases, we do not see evidence 

for resolved structure in H2 in any of the 23 targets. The velocity separation we could detect varies 
from target to target due to S/N issues and differences in strength and saturation level of the high- J 
lines. However, typical values arc 20-30 km s~^. 

We model the line-spread function of the spectrograph with a Gaussian of FWHM correspond- 
ing to a resolution, R ~ 17,000. This corresponds to the typical resolution of our spectra which were 
all observed through the largest available slit (30" x 30"). In many cases we are achieving greater 
resolving power, but even the J = 2 lines are considerably broadened beyond the instrumental 
profile, and the J = and 1 profiles are not affected by the choice of the line-spread function for 
any reasonable value. 

The H2 model itself includes the R(0), R(l), P(l), R(2), and P(2) lines of the band being fitted, 
as well as the R(0), R(l), and P(l) lines of adjacent bands. We must include the latter lines to 
account for the overlapping damping wings of the J = and 1 lines from adjacent bands. We use line 
parameters from Abgrall et al. (1993). Once this large model spectrum is calculated on a somewhat 
finer wavelength grid than the actual spectrum, it is convolved with the line-spread function, the 
zero-point shift is applied, and the model spectrum is rebinned to the grid of wavelengths from the 
actual spectrum. 

The final fits we report use the CURFIT routine with uniform weighting. In some cases the 
data quality is too poor at the shorter wavelengths and we can not adequately perform fits of the 
(4,0) band. Also, the SiC channels have poorer S/N in our range of interest, so in some cases we 
can obtain fits in the LiF channels but not SiC. A more subtle problem also occurs with the (4,0) 
band in certain cases. Considerable information is carried by the "bump" between the cores of 
the R(l) and P(l) lines. This small non-zero section in the spectrum at the saturated core of the 
vibrational bandheads is very sensitive to the J = 1 column density, and somewhat sensitive to 
the N{1)/N{0) ratio, yet it is usually weak enough to be insensitive to issues such as continuum 
placement. Due to the smaller /-values of the lines in the (2,0) and (1,0) bands, the bump is very 
prominent. But in the (4,0) band, a combination of large column density, poor data quality, and 
the larger /-values can totally eliminate the bump. This situation leaves the fitting routine with 
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only the blue and red wings of the overall J = and 1 profile to constrain the column densities 
and still accurately define the continuum. In some cases, this leads to unreasonable fits where the 
continuum will be strongly parabolic instead of relatively flat, and the column density may disagree 
with other bands by factors of 5-10. 

3.2. Error analysis 

In the ideal case we have 9 independent measurements of the column densities per observation. 
(For HD 73882 and HD 206267, the multiple observations give us more measurements.) Since the 
formal errors on the fits are smaller than the fit-to-fit differences, we adopt the same error analysis 
as in Paper II, and use as l-cr errors the sample standard deviation of the individual fits. For this 
choice to be appropriate, the individual fit parameters must be more-or-less normally distributed. 
We have 16 datasets for which we could obtain all 9 possible column density measurements, and we 
have used these datasets to search for systematic differences between the individual band/segment 
combinations. We determined the quantity 

^_ logiV(J),-(logiV(J)) 

for each of the measurements; i.e., the normalized deviation from the line-of-sight mean. If there 
are no systematic differences, the distribution of d should be normal. Then we found the average 
of the 16 d- values for each band/segment combination, along with the error of the mean. 

Table 7 gives the results of this analysis, along with the results for the same band in all 
segments, and all bands in all segments (whose average must be zero). We immediately see that 
there are a few systematic differences. The largest difference occurs for the LiF lA (2,0) .7 = fits 
which average nearly 1.5 standard deviations below the overall mean. The red wing of this profile is 
at the very edge of the detector segment which may affect the continuum determination. There are 
also relatively large systematic effects in the positive sense for both (1,0) band fits for J = 0, which 
effectively cancels out LiF lA (2,0) in the overall average. The source of these effects is consistent 
with the presence of a very weak stellar line in the blue wing of some of the (1,0) profiles, and an 
inspection of the fits does suggest that this is the case. In contrast with J = 0, the J = 1 values 
show much more subtle differences. 

Overall, it appears that while systematic differences do indeed occur, they do not have a 
large effect on the overall results. For example, if we were to exclude the LiF lA (2,0) fits from 
the averages, the logarithmic column densities would only change by a few hundredths. We also 
note that since we have excluded the (3,0) band we do not see band-to-band differences as large 
as those reported in Paper II. Interestingly, while the cases where we have the fewest individual 
measurements involve bands with large systematic differences (i.e., LiF lA (2,0), LiF lA (4,0), and 
LiF 2 A (1,0)), these systematic effects almost exactly cancel each other out when looking at the 
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ensemble for J = 0, and are small in any case for J = 1. For these reasons, we have included all of 
these bands in the final averages and uncertainties, i.e. the mean and sample standard deviation 
of the column densities from the (up to) 9 fits. 

4. Results 

Tabic 8 summarizes the measured and derived quantities relating to our H2 observations, which 
we will discuss individually in the following sections. We have generally used atomic hydrogen 
column densities from the literature, derived from lUE observations of the Lya line. In two cases, 
we report a new determination of iV(H I) from our own profile fits of lUE data. For stars later 
than about spectral type B2, the contribution from the stellar Lya line becomes large enough to 
seriously contaminate the interstellar line (§ 2.3.3). Thus, for six lines of sight we have estimated 
A^(H I) from the relationship between E{B - V) and A^(Htot) (see § 4.2). 

Our fundamental observed quantities are N{J=0) and N{J=1), from which we can derive 
the total molecular hydrogen column density (§ 4.1), total hydrogen column density (§ 4.2), kinetic 
temperature (§ 4.3), and hydrogen molecular fraction (§ 4.4). In addition, we can assess correlations 
between H2 parameters and extinction curve parameters (§ 4.5). The plots in these sections include 
lines of sight for which iV(H2) was measured by Copernicus (Savage et al. 1977). For these targets, 
we used N{}i I) from Diplas Sz Savage 1994 (lUE) and Bohlin et al. 1978 {Copernicus), in order of 
preference, with most values coming from the former. 

4.1. Molecular hydrogen column density 

While most of our targets have never been observed at moderate resolution in the far-UV, three 
targets in our present program were observed by Copernicus, providing measurements of A^(0) and 
A^(l), or at least A^(Il2). Table 9 compares those values with our new values. The difference for 
A''(l) for HD 24534 is quite large, but on the whole the differences are reasonable, given the much 
larger uncertainties on the Copernicus measurements due to poorer S/N. We also note that we have 
analyzed several Copernicus targets with our fitting techniques and find very close agreement with 
the published values. 

With the exception of the uncertain measurement of N(}l2) toward HD 24534 (X Per) by 

Mason et al. (1976), which we have refined, all but six of our present H2 column densities are 
larger than any observed with Copernicus. In four cases, our column density is larger than the 
revised value for X Per. We have thus provided the first significant sample of lines of sight with 
log iV(H2) ^ 21. 

Figure 2 shows A'^(Il2) as a function of color excess for log A'^(Il2) > 20. Although the column 
density appears to level off at large color excess, much of this leveling is due to our semi-logarithmic 



-23- 



presentation and the linear relationship between total hydrogen column density and color excess 
described in the following section. Interestingly, most of the scatter in A^(H2) for large color excess 
is due to scatter in A^(0), as A^(l) remains nearly constant at lO^*^'^ cm^^. For the 14 Ff/S'i? targets 
with E{B — V) > 0.4, the standard deviation of A^(0) is 0.21 dex as compared with 0.11 dex for 
A'"(l), and 0.15 dex for iV(H2). This finding is unchanged if we only consider values with uncertainty 
of 0.1 dex or less. However, it is not clear if this apparent threshold at 10^°-^ cm~^ is meaningful. 

Given our new H2 measurements, we can extend the range of column densities in exploring 
correlations with other molecules. In Figure 3 we show the relationships between N{}i2) and 
N^CR), A?^(CH+), A/"(CN), and N{CO). For the Copernicus H2 targets we took column densities 
from the compilations of Welty k Hobbs (2001) and Federman et al. (1994) for CH, Allen (1994) 
for CH"*", and Federman et al. (1994) again for CN and CO. We have only included absorption 
measurements since that step gives the best guarantee that we are sampling the same material as 
in the H2 measurements. 

The excellent relationship between CH and H2 seen in previous investigations continues to be 
reflected when adding our new data. The relationship is nearly linear, in agreement with Danks, 
Federman, & Lambert (1984). Chemical models predict a linear relationship between N{CH.) and 
A''(H2), with some scatter due to variations in density (Danks et al. 1984; van Dishoeck & Black 
1988, 1989). 

We also find a linear relationship between CH+ and H2, with some outlying points and consid- 
erable scatter. The primary formation reaction for CH"*" is endothermic, and thus shocks have been 
proposed as an energy source for the reaction (Elitzur & Watson 1978). Lambert & Danks (1986) 
found a good correlation between "warm" gas as measured by the rotational excitation of H2 for 
J=3-5, and A^(CH+), supporting the shock hypothesis. Our future measurements of rotationally 
excited H2 may shed additional light on this issue. 

Allen (1994) previously found that log A''(CH"'") increases with E{B — V) up to about 0.6, then 
levels off. However, this primarily occurs due to the semi-logarithmic axes and this relationship 
actually remains more or less linear up through E{B — V) = 1.2. Gredel (1997) also found linear 

relationships between A'^(CH"'") and extinction (Av) within individual OB associations. In addition, 
Gredel found a correlation between A^(CH+) and A^(CH). Gredel concluded that the dissipation of 
turbulence may be an important production mechanism. 

Another statistically significant relationship appears between CN and H2, similar to that found 
by Danks et al. (1984). The CN radical is highly density sensitive in the range of column densities 
studied here (Federman, Danks, & Lambert 1984). We also see a strong correlation between CN 
and molecular fraction (§ 4.2), and we would expect the latter quantity to also be correlated with 
density. The chemical models of van Dishoeck & Black (1988, 1989) appear to trace only the 
upper envelope of the observed CN abundances, but this model is for n ~ 500 cm~^. Low- and 
high-density models of Federman et al. (1984) appear to bracket our new data as they did for the 
Federman et al. dataset. 
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Finally, we see that CO and H2 are also highly correlated. The slope of the relationship appears 
to become more gradual at the highest column densities. At log A^(H2) ~ 20.5, the CO column 
density begins to increase more rapidly than H2. However, this is also the point where saturation 
effects become very important in assessing A^(CO) and in several cases a small 6-value (~1 km 
s~^) has been assumed which may not be appropriate if multiple components exist. Even a modest 
increase in the 6-value can result in a decrease in column density of an order of magnitude if the 
line lies on the fiat part of the curve of growth. 

On the other hand, chemical models do predict a rapid increase in CO column density relative 
to N{}i2) within the range of column densities studied in the present work. The models from 
van Dishoeck & Black (1988, 1989) show good agreement with the limited data. Absorption-line 

measurements of CO are difficult within the FUSE bandpass as the rotational structure is poorly 
resolved and the lines are often saturated. HST observations of the A-X series of CO lines will be 
crucial for extending the CO/H2 ratio. 

We note that while the correlations between the line-of-sight quantities are generally strong, 
some caution is necessary. As we discuss in § 5.2, we can not assess the true distribution of the 
majority of the H2 along the line of sight. Thus, in cases where the various molecules are distributed 
across a range of velocities, the line-of-sight correlations may not be physically meaningful. In 
particular, the column density of CN is sensitive to particle density and is expected to only trace 
the densest cloud cores. 

4.2. Total hydrogen column density 

First, we can look at the relationship between the total hydrogen column density, A^(Htot) = 
2A'"(H2) + A'^(H I), and color excess; i.e., the gas-to-dust ratio. Bohlin et al. (1978) found a linear 
relationship from Copernicus data, Af(Htot) = (5.8 x 10^-*^ cm~^) mag^^ E{B — V). In Figure 4, 
we have plotted the Copemicus/IUE dataset, along with our present FUSE sample. The new data 
fit the old relationship remarkably well, and the FUSE data alone give a slope of 5.6 x 10^^ cm~^ 
mag~^. There are a few disagreements larger than the error bars, but the largest deviation for 
E{B -V) > 0.3 is for the Copemicus/IUE observations of p Oph A {E{B -V) = 0.47, log (H I) 
= 21.7). This deviation still occurs with the revised lUEB. I measurement (Diplas &; Savage 1994) 
even though it represents a significant downward revision of the original Copernicus value. 

This enhanced gas-to-dust ratio for p Oph A has been interpreted as due to a preponderance 
of large grains within the p Oph cloud (Bohlin et al. 1978). These large grains are less efficient at 
producing visual reddening, and thus the E{B — V) color excess underestimates the actual quantity 
of dust. The unusual dust properties also lead to an unusual extinction curve and a very large value 
of Rv- We note that we have only a single line of sight in our present sample with Ry > 4 (HD 
102065) and we do not have an independent measurement of ^"(11 1) for this target. 
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4.3. Kinetic temperature 

In deriving the temperature of the gas, Tqi, we assume that the density and column density are 
high enough such that thermal proton collisions dominate over other processes in determining the 
ratio N{1)/N{0), and that the observed populations obey the Boltzmann relation. In textual form, 
we will refer to this temperature as the "kinetic temperature", while symbolically we will use Tqi 
to emphasize the source of this temperature. In addition, we emphasize that this temperature is 
a "column-averaged" temperature, while the actual temperature will vary throughout the cloud(s) 
in the line of sight. 

With a ratio of statistical weights, gi/go = 9, the population ratio is simply 



^(^) _ Op-Eoi/feToi /o^ 
N{0) ~ ' ^ ^ 

where Eoi/k = 171 K. With column densities expressed as base-10 logarithms (as in Table 8), the 
kinetic temperature (in K) can then be written 

74 

" log iV(0) - log Ar(l) + 0.954' 

In calculating the uncertainties in kinetic temperature, we take the combination of 1-a errors 
that gives the largest deviation from the best value, such that the errors on the derived values are 
more conservative. Furthermore, in deriving these errors, we have taken 0.04 dcx as the minimum 
possible error on a column density, even when we have derived a smaller error. This corresponds 
to 10% and while this choice is arbitrary we feel that it is a reasonable guess for the magnitude of 
any systematic effects. 

The average kinetic temperature derived from Copernicus observations of 61 lines of sight with 
log A''(H2) > 18.0 was 77 ± 17 K (Savage et al. 1977). A similar calculation for the 9 Copernicus 
lines of sight with log A^(H2) > 20.4, comparable to the present survey, gives 55 it 8 K. Our FUSE 
sample gives an intermediate value, 68 it 15 K. However, we note that our sample has a somewhat 
unusual distribution, with three lines of sight having Tqi > 94 K, but none in the range 75-93 K. 
In any case, our average value is similar to that found previously for lines of sight where H2 is 
self-shielded. 

Despite extending the range of color excess by a factor of 2, Figure 5 shows that the kinetic 
temperature in our sample does not change with increasing E{B — V). We have also searched for 
a correlation between Tqi and Ry and found none. We might expect the kinetic temperature to 
be anti-correlated with density indicators, and we see such a relationship between Tqi and A'^(CN) 
(Figure 6). The slope of the relationship is quite small and there are a few outlying points, but 
given the small range in the observed temperatures, the relationship is quite good. 
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4.4. 



Molecular fraction 



The hydrogen molecular fraction, /h2, gives the fraction of hydrogen atoms in molecular form. 
In terms of the column densities of H I and H2 , 



The tabulated uncertainties for molecular fraction follow the same procedure described in the 
previous section. 

The Copernicus data showed an interesting trend of molecular fraction with increasing color 
excess (Savage et al. 1977). Below E{B — V) 0.08, the molecular fraction is quite small, typically 

less than 10^^, while above this point the fraction is generally greater than 10~^, with few points 
lying in between. This abrupt boundary occurs due to increased self-shielding of H2 near iV(H2) ~ 
10^^ cm"^, corresponding to /h2 ~ 10~^. 

In Figure 7, we show molecular fraction versus color excess for the Copernicus data and our 
FUSE data. The boundary at E{B — V) ^ 0.08 is not visible because we have chosen a linear 
scale for the ordinate. The FUSE data mostly overlap values found previously, but we have greatly 
increased the number of lines of sight with at least moderately high molecular fraction. Even in 
the range of overlap of the two samples, the FUSE sample shows larger molecular fractions, but 
this is probably a selection effect. We do not see an increase in molecular fraction with increasing 
extinction within the FUSE sample. 

Figure 8 shows the molecular fraction versus the total-to-selective extinction ratio, Ry, for the 
FUSE dataset. Previous results have suggested an anti-correlation between the two quantities for 
diffuse clouds, consistent with idea that grain coagulation reduces the available surface area for H2 
formation at larger Ry (Cardelli 1988). However, our data do not show a statistically significant 
relationship between the two quantities. We do not presently have good coverage of large values of 
Rv, but several lines of sight with Ry > 4 remain to be observed as part of the FUSE translucent 
cloud program. 

Although little-mentioned, the Copernicus data show a good correlation between molecular 
fraction and kinetic temperature. Figure 9 shows both data sets, and while the FUSE data them- 
selves show only a weak relationship, those points still follow the general trend. There are lines of 
sight with small molecular fraction at all kinetic temperatures, but lines of sight with large molec- 
ular fraction are preferentially associated with small kinetic temperature. This relationship is not 
surprising as both large molecular fraction and small kinetic temperature should be associated with 
denser cloud cores. In a similar sense, we also see a correlation between molecular fraction and the 
density-sensitive CN abundance (Figure 10). 



/h2 = 



2N{H2) 



(10) 



2N{H2) + N{HI) 
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4.5. Extinction curve parameters 

With six extinction curve parameters and four column-density related quantities {fti2, -^(CH)/A^(Htot)) 
Af(CH+)/Af(Htot), A^(CN)/iV(Htot)), we have 24 potential correlations. In addition to our FUSE 
data points, we have included the handful of points from the Fitzpatrick & Massa (1986, 1988, 
1990) and Jcnniskcns & Grccnburg (1993) extinction curve surveys for which we have the ancillary 
data. We have evaluated the Spearman rank correlation coefficient for each of these relationships. 
In 11 cases the correlation is not significant at the la level, while 10 correlations are significant at 
the 2(7 level. Thus, we generally either find no correlation or a good correlation in the statistical 
sense. The 2a group includes N (CN) / N (Htot)) vs. both a and C2; N{CE.+)/ N{Rtot) vs. C3; fm, 
A^(CH)/A^(Htot), and N (CN) / N (Rtot)) vs. both C4 and 7; and A^(CH)/A^(Htot) vs. Aq ^ Since the 
main focus of this paper is H2, we will consider the two strong correlations involving /h2 in detail. 

The strongest {3.7a) and most intriguing correlation is that between molecular fraction and 
the width of the 2175 A bump, 7 (Figure 11). In fact, of all the parameters we have considered in 
the present work, 7 appears to be the best predictor of molecular fraction. The larger molecular 
fractions appear to be associated with regions of larger density. Thus, our findings suggest that the 
width of the 2175 A bump is closely related to density. Fitzpatrick & Massa (1986) reported a similar 
finding in a qualitative sense; dense quiescent regions such as dark clouds and reflection nebulae 
were associated with broad bumps, and diffuse clouds and star-forming regions were associated 
with narrower bumps. They found a good correlation between 7 and E{B — V)/r, where r is the 
distance to the star, even with the biases associated with this density indicator. Our observed 
correlation between /h2 and 7 shows much less scatter and could well be taken as linear. The 
single outlying point, ^ Oph, showed the largest difference in 7 between the two methods used by 
Fitzpatrick & Massa (1986, 1990). We have used the final value of Fitzpatrick & Massa (1990) 
based on the overall extinction curve fits, but the initial value, based on just the region around the 
bump, would lie much closer to the rest of the points. The authors attributed this large difference 
to the relative shallowness of the bump and the resultant uncertainty in separating the bump from 
the rest of the extinction curve. 

Ignoring ( Oph, an unweighted linear fit of the rest of the points gives the relation, 

/h2 = 1.447 - 1-04 (11) 

corresponding to a minimum value of 7 of 0.72, and a maximum value of 1.42. (If we include 
Oph in the fit, the allowed range is 0.67-1.51.) In fact, the extrema in 7 in the entire Fitzpatrick 
&; Massa (1990) and Jenniskens & Greenberg (1993) samples of more than 100 curves are 0.76 and 
1.383 (or 1.25 if we ignore C Oph). 

While the bump width appears to be a good predictor of molecular fraction in our sample, 
this relationship may not hold in all environments. For example, many lines of sight in the SMC 
have no discernible 2175 A bump at all (Gordon &; Clayton 1998). On the other hand, in a survey 
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of 30 Galactic lines of sight selected to sample low-density gas, Clayton, Gordon, & Wolff (2000) 
found small values of 7 (^1.0). This finding is consistent with the likelihood that the lines of sight 
sample gas with low molecular content. Thus, a strong correlation between molecular fraction and 
bump width may apply in most Galactic environments, albeit the relationship may not be linear. 

The strength of the far-UV curvature, C4, exhibits the other strong correlation with /h2, at 
the 2.6a level (Figure 12). The extrema in C4 are HD 102065 and HD 62542, and these points 
show up as outliers in the relationship. Although neither line of sight has a well-determined value 
of /h2 it is highly unlikely that our reported values are so inaccurate as to match the apparently 
linear trend with C4 seen in the other points. The presence of this correlation is consistent with the 
previously noted tendency for a steep far-UV rise to be associated with a broad 2175 A and the 
observed correlation between the bump width and molecular fraction. 

An alternative explanation for the correlations between molecular fraction and both the bump 
width and far-UV curvature concerns the properties of the dust grains. Increased far-UV curvature 
is thought to be associated with smaller than normal dust grains (Cardelli, Clayton, &; Math- 
ias 1989). The 2175 A bump is most likely associated with small carbonaceous grains (Desert, 
Boulanger, & Puget 1990), and perhaps smaller grains lead to broader bumps. 

Grain size is also thought to be smaller in lines of sight with small Ry, and Cardelli (1988) 
found an inverse correlation between molecular abundances and Ry- He attributed this correlation 
to the effects of these smaller grains and their effect on H2 formation and destruction. With 
similar total grain masses, the smaller grains will provide greater surface area, yielding a greater 
H2 formation rate, and a smaller photodissociation rate via the increased far-UV extinction. We 
recall, however, that in the present work, we do not find a good correlation between molecular 
fraction and Ry- 

5. Discussion 

In this section, we will mainly focus on what our findings say about the nature of our present 
lines of sight relative to diffuse clouds. The overall line-of-sight characteristics of most of our present 
sample satisfy the criterion to be considered "translucent", i.e. Ay > 1. Implicit in the definition 
of a translucent cloud is that we are considering a single molecular cloud, and not a collection of 
several diffuse clouds. For the purposes of this discussion, we adopt a definition of a "translucent 
cloud" similar to that envisioned by van Dishoeck &; Black (1988); i.e. /h2 ^ 0.9, Tqi < 40 K, and 
Ay > 1. Such a cloud may be an isolated cloud, a skin around a dense cloud, or a core located 
within significant diffuse material. 

If a line of sight is dominated by one of these clouds, we would expect this situation to be 
reflected in several of our measured quantities. Specifically, the observed molecular fraction should 
be large, while the kinetic temperature should be small. As shown in Figure 9, these two quantities 
do indeed show an anti-correlation, with considerable scatter. Despite the scatter, all of the lines of 
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sight where N{'K2) > N{]i I) (/h2 > 2/3) show small kinetic temperatures. In fact, there appears 

to be a distinct group of 10 lines of sight centered near Tqi = 55 K and /h2=0.7 that is separated 
from the rest of the samplc^^. This is even more apparent if we ignore the FUSE data points 
without direct measurements of N(H I). 

However, there are several lines of sight with similar or even smaller kinetic temperatures than 
this group of 10. In addition, while the molecular fractions are relatively large, they do not closely 
approach unity as we might expect. Also, the extinctions for several of these lines of sight are less 
than or equal to one magnitude, barely satisfying the rather loose definition of a translucent cloud 
we have adopted. To further assess the question of whether we are seeing individual translucent 
clouds, we need to consider the distribution of material along the hue of sight (§ 5.1). We also 
consider evidence from studies of chemical depletions (§ 5.2). Finally, we consider the question of 
why we see few, if any, translucent clouds in our lines of sight (§ 5.3). 

5.1. Multiple clouds and "hidden" translucent clouds 

Arguing in favor of the hypothesis that we are seeing at least a few translucent clouds is 
the fact that the overall line-of-sight column densities can be greatly affected by the particular 
distribution of material. Even if highly molecular material exists, there could be a skin of diffuse 
material surrounding this cloud, or additional diffuse clouds along the line of sight. In these cases, 
the observed integrated molecular fraction could be considerably less than unity and the kinetic 
temperature could be affected as well. Thus, even a line of sight with only a moderately high 
molecular fraction could harbor a translucent cloud. 

Chemical modeling of the lines of sight can, in principle, help constrain the distribution of 
material. Models have been constructed which reproduce the column densities of a variety of 
species in diffuse clouds, including H2 (e.g., van Dishoeck & Black 1986). In some lines of sight, 
the models have difficulty reproducing the high- J column densities of H2 (e.g.. Paper II). 

Part of our group has developed a new code to study the formation, destruction, radiative 
transfer, and ro- vibrational excitation of H2 (Browning et al. 2002). Several lines of sight with 
published high- J column densities for H2 have been modeled with this code, including HD 110432 
and HD 73882. In many cases, it was difficult to model a line of sight as a single cloud without 
invoking extreme conditions, such as a radiation field 10-100 times the Galactic mean, which is 
unrealistic in most cases. Better matches require multiple clouds and/or multiple pathways for 
incoming UV radiation. In addition, changes in one physical parameter can mask changes in 
another such that a group of models corresponding to a range of physical conditions will all match 
the column densities for an individual line of sight. Ensembles of models, combined with large 



i^These 10 lines of sight arc: HD 24534, HD 27778, HD 62542, HD 73882, HD 99675, HD 154368, HD 210121, and 
the Copernicus targets C Oph, o Per, and C Per 
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samples of lines of sight, can provide more definitive results, such as in the LMC/SMC (Tumlinson 
et al. 2002). 

We can also observationally study the possibility that the H2 along these lines of sight lies 
in several distinct clouds, with corresponding smaller extinctions for each cloud. The very high 
resolution ground-based data show that most lines of sight have numerous components in Na I and 
K I, and the strongest components usually lie within a 10 km s~^ velocity range (D. E. Welty et al. 
2002, in preparation). Data for the CH radical show similar behavior. Given the good correlations 
between N{}i2) and these species, we might expect the H2 to also be distributed among several 
components. 

We can assess the distribution of H2 indirectly, in several ways. Equivalent width data allow 
us to perform a curve-of-growth analysis of the high-J lines. Most simply, we can construct a 

best-fit single-component curve of growth via fitting techniques and determine the best values for 
the column densities of each J-lcvcl as well as the optimum "effective" 6-valuc of the distribution. 
However, we can not always assume that the actual component structure is well approximated by 
a single 6-value. Our program of high-resolution spectroscopy allows us to better constrain the 
component structure with the assumption that some or all of the components seen in Na I and/or 
K I and/or CH contain H2, with perhaps a scaling of the 6- values of the individual components, and 
possible variation in the ratios of the ground-based species to H2 from component to component. 
As we have seen, for the overall lines of sight, there is a nearly one-to-one correspondence between 
CH and H2, and strong correlations also appear for Na I and K I (Welty & Hobbs 2001). 

In our preliminary investigations of the high-J H2 lines, we often find a better match with 
"effective" 6-values greater than can be obtained with the CH component structure. We would 
then have to use additional components observed in K I or Na I — including some of the weaker 
components, as compared with only using the strongest. Thus, even with a "simple" component 
structure, the H2 still appears to be distributed among several closely-spaced velocity components. 
Despite the small separations, the components arc resolved in the very high resolution optical spec- 
tra, which suggests discrete components and not just material more-or-less uniformly distributed in 
velocity space. With the additional components that may contribute at least small amounts of H2 
to the overall observed column densities, the quantity of H2 available in each component is further 
diluted. 

One important caveat to this discussion of the component structure is that in looking at the 
high-J lines we are considering only 1% or less of the total H2. We have little hope of directly 
assessing the true component structure of the extremely strong J = 0-1 lines. In addition, the J 
= 2 lines are also very strong and heavily saturated or even damped, and will not always provide 
constraints to the component structure. We expect the material containing the low- J lines to mostly 
be located in the self-shielded cloTid core(s), while the material containing the high-J lines may be 
more physically widespread along the line of sight. Thus, the component structure corresponding 
to J = and 1 may be simpler than that for the high-J lines. 
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There is a tendency for the Hnes of sight with the largest molecular fractions to have the 
simplest component structures when looking at CH or the strong Na I and K I components that 
are most likely to also contain H2 (D. E. Wclty et al. 2002, in preparation). In Table 10 we give 
preliminary information on the CH component structures observed at high resolution. We caution 
that the 4 cases with a single component are for observations at the lowest resolutions. The typical 
spacing between components is 3-5 km s~^, so there may be unresolved structure, as indicated 
in higher resolution K I spectra in these 4 cases. In fact, we find more complex structure in the 
K I spectra in all lines of sight in Tabic 10. It is more difficult to detect weak components or 
resolve closely spaced components in CH observations due to A-doubling, the small atomic weight, 
and the relatively weak spectral lines. In the CH data with resolution of less than 2 km s~^, 
the lines of sight with more material in the strongest component and/or fewer components have 
larger molecular fractions. These lines of sight presumably have less contamination from foreground 
clouds of atomic gas. Yet, in all cases where we have this "very high" resolution data, there are 
multiple components, with no more than ^-^75% of the molecules in the dominant component. 

Given the relatively close spacing of the CH components we have to consider whether they 
represent small knots of material within a single cloud, or distinct clouds along the line of sight. 
If we are dealing with knots of material within a single cloud, we still might consider this cloud 
to be a translucent cloud due to the large overall extinction and high molecular content. Also, 
within a single cloud, the extreme widths of the H2 J = and 1 lines effectively self-shield the 
entire range of velocities. If we are dealing with multiple distinct clouds at the even closer spacing 
suggested by the K I data, the amount of molecular material within the most abundant cloud core 
will generally be even less than indicated in Table 10 (or in our mathematical analysis below and in 
Table 11). This suggests that these lines of sight are simply collections of many diffuse clouds with 
small extinction, albeit with large molecular fractions and presumably relatively high densities. 

An indirect method for assessing the quantity of H2 and H I in each component comes from 
the study of chlorine. The chemistry of chlorine is intimately connected to that of H2 because 
of the large reaction rate with CI II. With a first ionization potential of 13.01 eV, chlorine will 
be primarily ionized in H I gas, and primarily neutral in association with H2 (Jura 1974; Jura k. 
York 1978). Thus, observations of CI I and CI II lines (in the UV) can constrain the hydrogen 
molecular fraction. P. Sonnctruckcr ct al. (2002, in preparation) have used this technique to infer 
that the true molecular fraction in the dominant area of H2 toward HD 192639 is as large as twice 
the line-of-sight value (2/3 vs. 1/3). High-resolution measurements of the the CI 1 A1347 line may 
provide important constraints on the molecular fractions of individual velocity components toward 
our lines of sight. 

While detailed physical and chemical modeling of these lines of sight will appear in a later 
paper or papers, we can perform a feasibility analysis on the H2 column densities in each line of 
sight to assess the possible presence of translucent clouds. I.e., we know the overall line of sight 

parameters, and we also know approximately what a translucent cloud should look like. Thus, we 
can add the expected observed properties of putative translucent clouds to diffuse cloud material in 
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various combinations that produce the observed total hne of sight quantities. While not a physical 
model of the clouds, it does provide interesting constraints on the observable quantities. 

We define translucent material as having /h2 = 0.9 and Tqi = 30 K, comparable to the param- 
eters used in the translucent cloud models of van Dishoeck & Black (1988). To constrain the total 
column density of the translucent cloud and the remaining material, we assume that the previously 
discussed relationship between E(B — V) and Af(Htot) still applies for highly molecular clouds. We 
then mathematically create the largest embedded translucent cloud that is still consistent with the 
line of sight totals, and that produces realistic parameters for the remaining "diffuse" material, i.e. 
Toi < 200 K. If the line of sight molecular fraction is small and the temperature is large, only a 
small amount of material could possibly be locked up in a cold, highly molecular cloud, and vice 
versa. 

In Table 11 we give our results for this breakdown of material into translucent and diffuse 
clouds. From this simulation, we see little evidence for individual translucent clouds along our lines 
of sight. In two directions, HD 73882 and HD 154368, translucent clouds with Ay ~ 1 could exist. 
The remaining diffuse material still has a molecular fraction appropriate for the color excess. HD 
62542 is the only direction where more than half of the material could be translucent. 

If we were to relax the criteria for a "translucent" cloud to something like /h2 = 0.8 and Tqi 
= 40 K, comparable to the most extreme line of sight values we actually see, we would then have 
many more cases where there is as much translucent material as diffuse. Not coincidentally, this 
larger list of lines of sight with as least as much translucent material as diffuse is identical to the 
list of FUSE targets in the upper left-hand portion of Figure 9 discussed at the beginning of § 5. 
However, we would still not have additional translucent clouds with Ay > 1. 

We have assumed a gas-to-dust ratio of 5.8 x 10^^ cm mag~^ based on the excellent match 
we see in Figure 4. The observed scatter could support variations on cither side of this relationship 
for our FUSE targets, albeit not as large as seen for p Oph. However, we would need large (and 
unrealistic) differences in the dust characteristics between the "diffuse" and "translucent" material 
to modify our conclusions in the previous paragraphs. 

5.2. Evidence from chemical depletions 

As noted in § 2.3.13, chemical depletions can be used as an indicator of the presence of translu- 
cent clouds. For HD 192639, the depletions of about a dozen species indicate similar conditions 
to diffuse clouds (P. Sonnentrucker et al. 2002, in preparation). We have also undertaken a study 
of Fe II depletions in most of our FUSE targets (Snow, Rachford, & Figoski 2002). Snow et al. 
found that the depletion of iron is more or less uniform for the FUSE lines of sight. There is 
little evidence for increased depletion with increasing extinction, molecular fraction, or A?^(Htot)/^, 
within the range of E{B — V) and Ay covered by FUSE. The new observations did not extend 
the trends of increasing depletion with increasing density found by Savage &; Bohlin (1979) and 
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Jenkins, Savage, & Spitzer (1986). Such increased depletions are expected to occur at some point 
within the translucent clouds (Snow et al. 1998). Given the moderate resolution of the FUSE 
Fe II observations, wc can not rule out the possibility of large depletions for individual components, 
which become "smoothed out" in the line-of-sight values. 

Two of our Ff/S'-E targets with particularly large molecular fraction, HD 24534 and HD 154368, 
have been observed at high resolution with HST to search for increased depletions in several addi- 
tional species. The results for HD 24534 show a very slight, and probably not statistically significant, 
increase in carbon depletion relative to diffuse clouds (Sofia, Fitzpatrick, & Meyer 1998). Snow et 
al. (1998) did not find increased oxygen depletion toward HD 24534. Snow et al. (1996) studied 
many species toward HD 154368 and also did not find increased depletions. They concluded that 
the line of sight toward HD 154368 contained extended regions of moderate-density gas instead of 
one or more dense cloud cores, despite Ay ~ 2.5. 

5.3. Where are the translucent clouds? 

While each of techniques we have discussed has limitations, we are forced to conclude that with 
few exceptions there is little evidence for individual translucent clouds within our lines of sight, 
based on the definition given at the beginning of this section. For HD 24534, HD 154368, and HD 
192639 we have considerable information that leads to this conclusion, while the evidence for the 
rest of the lines of sight is more circumstantial. There are several possible reasons why we have not 
found such clouds: 

1. Translucent clouds are just beyond our FUSE range. 

Very few stars have enough UV flux to be observed with FUSE through Ay > 2, and our present 
sample includes only seven such targets. Many interesting targets were rejected from the original 
potential target list for being too faint. Several of these targets have Ay = 2-4 (e.g. HD 80077, 
HD 169454) and lie behind molecular clouds. When combined with the fact that many, if not most, 
lines of sight will have contaminating diffuse material, we may need to reach Ay ~ 3-5 to find 
evidence for individual clouds with molecular fractions near unity and Ay ~ 1-2. 

Chemical models suggest that several important transitions involving carbon take place in 
clouds in the range Ay ~ 1-3, such as increased CO abundance and a dramatically larger C I/C II 
ratio (e.g., van Dishoeck & Black 1988, 1989). There is also evidence for changes in grain conditions 
at Ay ~ 3. For example, Whittet et al. (2001) find an increase in Ry in the Taurus dark clouds 
for Ay > 3. The 3-/xm water-ice feature begins to appear at this point as well, lending support for 
the idea that grain mantle growth becomes important here. 

2. ^ Oph-type lines of sight are indeed "translucent." 

As we have mentioned, large molecular fractions do seem to be correlated with not only simpler 
component structures, but also low kinetic temperatures, large column densities of density-sensitive 
species such as CN, and the environment-sensitive 2175 A bump width and far-UV extinction 
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curvature. If we were to relax the requirement for large Ay, we could consider all 10 lines of sight 

discussed in the preamble to this section as containing "translucent clouds." However, we note 
that even within this group, the already small Ay is sometimes broken up into even smaller clouds 
(such as toward C, Oph) . We are reluctant to accept this possibility without a confirmation that we 
do not see large molecular fractions in appropriate targets with large Ay. 

3. There is no distinction between "diffuse" and "translucent" except larger Ay. 
This possibility is the simplest. However, there is significant evidence listed above to support the 
idea of at least a gradual transition in physical parameters through Ay = 1-5, if not more abrupt 
transitions. Again, we feel that we need to cover a larger range of Ay before we could accept this 
possibility. 

We feel that the first possibility is the most likely one based on our present findings. The 

ubiquity of interstellar clouds makes it difficult to find lines of sight that only sample a single 
cloud, or stars that are still observable through an ensemble of several clouds that might contain a 
translucent cloud. 

A more definitive answer to the question posed in this section would require H2 observations of 
several lines of sight with Ay > 2 where the neutral atoms and simple carbon-containing molecules 
show relatively simple velocity structure such as HD 169454 (Jannuzi et al. 1988; Crawford 1997). 
Even in this case, the Na I component structure is more complicated (Federman Sz Lambert 1992). 
A more serious problem is that FUSE is barely capable of this observation, and is incapable of 
observing several other promising targets such as HD 80077. We note that while the decrease in 
observed UV fiux due to increasing H2 absorption plays a role in limiting the number of targets 
we can access, it is a very small role. For instance, we can calculate model H2 absorption spectra 
and compare the height of the peaks between the vibrational bandheads for N(H2) = 1 x 10^^ 
(present sample) and 5 x 10^^ cm~^ (a highly molecular cloud with Ay 5-6). For this extreme 
example, the intensity relative to the continuum at the peaks on either side of the (4,0) bandhead 
is ~0.85 for the low-column case, and ~0.45 for the high-column cases. However, the effect at the 
(2,0) and (1,0) bandheads is much smaller since these bands arc weaker, i.e., ~0.95 versus ~0.75. 
This is a trivial effect relative to the increased dust extinction, and indicates that our dataset is 
not significantly biased against high-column lines of sight through selection effects. 

As a final note, we emphasize that there is a subset of lines of sight with relatively small 
extinction, but extreme line of sight characteristics. These include HD 62542 and HD 210121, both 
of which appear to have a large molecular fraction, and are among the best candidates for having 
translucent clouds based on Table 10. We may eventually find lines of sight with larger Ay that 
indeed contain individual clouds similar to those toward HD 62542 and HD 210121 and thus the 
latter would be considered translucent clouds. 
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6. Summary 

We have completed the first FUSE survey of molecular hydrogen in lines of sight with Ay ^ 
1. The survey includes observations toward 23 early- type stars, which sample gas in a range of 
environments. Through profile fitting of three vibrational bandheads, we have directly measured 
the H2 column densities in the J = and J = 1 states. Combined with ancillary data, we have 
then derived the total H2 column density, hydrogen molecular fraction, and kinetic temperature. 
In addition, we have compiled a set of extinction parameters for our lines of sight. 

With this information we have investigated many important relationships between parame- 
ters. We have extended previous correlations between A^(H2) and A^(CH), A^(CH+), A^(CN), and 
A^(CO). While the formation of CH+ is still an open question, the other correlations closely match 
predictions from chemical models, despite possible differences in the distribution of these species. 
We find a potentially powerful combination in assessing the H2 content of a line of sight based 
on ground-based measurements and lUE or HST mid-UV observations. The excellent relationship 
between A^(H2) and A^(CH) gives a good estimate of A^(H2), while the measurement of the width 
of the 2175 A bump gives a good estimate of molecular fraction. 

We find a generally self-consistent picture of these lines of sight in the sense that various 
indicators of density correlate well with each other; i.e., molecular fraction, kinetic temperature, 

CN abundance, and extinction curve parameters. 

While our sample has a relatively large average molecular fraction, we have not found lines of 
sight with molecular fractions greater than 0.8. We have identified a subset of 7 lines of sight (plus 
3 lines of sight observed by Copernicus) with large molecular fraction and small kinetic temperature 
which represent the best candidates for the presence of "translucent cloud" material. However, the 
possible quantities of material that could be associated with translucent clouds generally corre- 
sponds to Ay < 1 mag, the nominal lower limit for such clouds, and in no case do we see evidence 
for highly molecular material corresponding to Ay ^1-5 mag. In addition, in most cases we see 
evidence for multiple velocity components for H2, which would further divide up the extinction 
when considering the individual clouds. Thus, our conclusion is that for the present sample we are 
observing only a few individual translucent clouds. Rather, we are mostly seeing combinations of 
diffuse clouds. We suggest that without specific evidence to indicate individual translucent clouds, 
lines of sight with Ay > 1 should be called "translucent lines of sight" and the term "translucent 
cloud" should be avoided. 

An ongoing analysis of the high- J lines of H2 in these lines of sight, combined with a detailed 
modeling program similar to that for the Magellanic Clouds (Tunilinson ct al. 2002), will allow 
us to better assess the physical conditions in the clouds. This will be the first large study of 
high-J excitation of H2 to incorporate very high resolution ground-based observations of species 
such as K I, Na I, and CH to help assess the unresolved component structure of H2. Also, we are 
undertaking detailed analysis of all available species in selected individual lines of sight (e.g. P. 
Sonnentrucker, et al. 2002, in preparation, for HD 192639), combining FUSE data with recently 
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obtained HST data. Further studies with HST in the hues of sight with the largest molecular 
fractions will also be useful. 

Finally, several additional lines of sight from the FUSE translucent cloud survey have already 
been observed since the cut-off date for inclusion in the present work (June 2001), and we anticipate 
that most of the 21 remaining targets will be observed. Several very important targets such as HD 
37903 and HD 147889 are part of the program. These additional observations will double the 
number of lines of sight studied, and should improve and/or extend the results of the present work. 

We thank the anonymous referee for many helpful comments. This work is based on data 
obtained for the Guaranteed Time Team by the NASA-CNES-CSA FUSE mission operated by the 
Johns Hopkins University. Financial support to U.S. participants has been provided by NASA 
contract NAS5-32985. This research has made use of the SIMBAD database, operated at CDS, 
Strasbourg, Prance. 
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1050 1060 1070 1080 

Wavelength (A) 

Fig. 1. — Sample spectra at high and low S/N. Top - HD 210839 (peak S/N 24 per pixel); 
bottom - HD 154368 (S/N ~ 2). Both spectra have been smoothed to resolution (~9 pixels) for 
presentation. This spectral region includes the Lyman series (4,0), (3,0), and (2,0) vibrational 
bandheads of H2 (from left to right), as well as broad stellar features at 1057 A and 1073 A. A less 
obvious stellar feature lies near 1067 A. The numerous narrow lines are mostly due to rotationally 
excited H2, as indicated by the series of tickmarks above the spectra; exceptions include Ar I A A 
1048,1067 and Fe II A1055. Profile fits to the (4,0) and (2,0) bands are shown. 
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Fig. 2. — Molecular hydrogen column density vs. color excess. Symbols - crosses: FUSE; diamonds: 
Copernicus 
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Fig. 3. — H2 column density vs. other molecular column densities. Symbols as in Figure 2. Solid 
curves in panels a), b), and d) are translucent cloud models from van Dishoeck & Black (1989). 
Dotted lines in the bottom right corner of each panel correspond to unit slope for that panel. 




Fig. 4. — Symbols as in Fig. 2. The six FUSE points with no independent measurement of 
N{}i I) are not included. The solid line corresponds to the relation A?^(Htot) = (5.8 x 10^^ cm~^ 
mag-^)E{B - V) given by Bohlin et al. (1978). 
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Fig. 5. — Kinetic temperature vs. color excess. Symbols - crosses: FUSE; diamonds: Copernicus 
points with iV(H2) > 10^° cm ^; squares: Copernicus points with A^(H2) < 10^° cm ^. 
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Fig. 6. — Kinetic temperature vs. fractional CN abundance. Symbols as in Figure 5. (Note the 
particularly large error bar for BD +31° 643 on this expanded vertical scale.) 
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Fig. 7. — Molecular fraction vs. color excess. Symbols - crosses: FUSE; asterisks: FUSE points 
with no independent measurement of N(H I); diamonds: Copernicus points with N{'H.2) > 10^° 
cm~^; squares: Copernicus points with N{}i2) < 10^° cm~^. 
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Molecular fraction vs. total-to-selective extinction. Symbols as in Figure 5. 
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Fig. 9.— 



Molecular fraction vs. kinetic temperature. Symbols as in Figure 7. 
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Fig. 10.— 



Molecular fraction vs. fractional CN abundance. Symbols as in Figure 7. 
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Fig. 11. — Molecular fraction vs. 2175 A bump width. Symbols as in Figure 7; note though that 
there are no Copernicus points with N(H.2) < 10^^ cm~^. 
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Fig. 12. — Molecular fraction vs. far-UV extinction curvature. Symbols as in Figure 7; note though 
that there are no Copernicus points with A'^(H2) < lO^'^ cm~^. 
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Table 1. Target list 



Star 


e 




b 




Assoc. 


V 


MK type 


Ref 


BD ^ 


h31° 643 


160, 


,49 


-17. 


,80 


Per 0B2 


8.51 


B5 V 


1 


HD 


24534 


163, 


.08 


-17. 


,14 


Per 0B2 


6.10 


BO Vc 


2 


HD 


27778 


172, 


.76 


-17. 


,39 


Tau-Aur 


6.49 


B3 V 


3 


HD 


62542 


255, 


.92 


-9. 


,24 




8.04 


B3 V 


4 


HD 


73882 


260, 


.18 


+0. 


.64 


Vel OBI 


7.24 


08.5 Vn 


5 


HD 


96675 


296, 


.62 


-14. 


,57 




7.67 


B7 V 


6 


HD 


102065 


300, 


.03 


-18. 


,00 




6.61 


B9 IV 


7 


HD 


108927 


301, 


.92 


-15. 


,36 




7.78 


B5 V 


7 


HD 


110432 


301, 


.96 


-0. 


,20 




5.24 


Bl Ille 


8 


HD 


154368 


349, 


.97 


+3. 


,22 




6.13 


09.5 lab 


5 


HD 


167971 


18, 


.25 


+1. 


,68 


Ser 0B2 


7.50 


08 Ib(f)p 


9 


HD 


168076 


16, 


.94 


+0. 


,84 


Ser OBI 


8.21 


04V((f)) 


5 


HD 


170740 


21, 


.06 


-0. 


,53 




5.76 


B2 V 


10 


HD 


185418 


53 


.60 


-2. 


,17 




7.45 


B0.5 V 


11 


HD 


192639 


74, 


.90 


+1. 


,48 


Cyg OBI 


7.11 


07Ib(f) 


12 


HD 


197512 


87, 


.89 


+4. 


.63 


Cyg 0B7 


8.57 


Bl V 


12 


HD 


199579 


85, 


.70 


-0. 


,30 


Cyg 0B7 


5.96 


06 V((f)) 


12 


HD 


203938 


90, 


.56 


-2. 


,33 


Cyg 0B7 


7.09 


B0.5 V 


12 


HD 


206267 


99 


.29 


+3, 


,74 


Cep 0B2 


5.62 


06.5V((f)) 


12 


HD 


207198 


103, 


.14 


+6. 


,99 


Cep 0B2 


5.96 


09.5Ib-II 


12 


HD 


207538 


101, 


.60 


+4. 


,67 


Cep 0B2 


7.31 


09 V 


12 


HD 


210121 


56, 


.88 


-44. 


,46 




7.67 


B3 V 


13 


HD 


210839 


103, 


.83 


+2. 


,61 


Cep 0B2 


5.05 


06 Infp 


12 



References. — (1) Harris, Morgan, &: Roman 1954; (2) Roche et al. 
1997 (3) Kenyon et al. 1994; (4) Whittet et al. 1993; (5) Walborn 1973; 
(6) Vrba & Rydgren 1984; (7) Houk 1975; (8) Slettebak 1982; (9) Walborn 
1972; (10) Murphy 1969; (11) Abt 1985; (12) Garmany & Stencel 1992; 
(13) Welty & Fowler 1992. 



-55- 



Table 2. Carbon-based molecular abundances 



Star 


log N{CR) 
(cm-2) 


Ref 


log iV(CH+) 
(cm-2) 


Ref 


log Ar(CN) 
(cin-2) 


Ref 


log Ar(CO) 
(cm-2) 


Ref 




16. 1 1 


i 


1 Q no 


i 


<iz.U4 


-1 
i 






xlU z4oo4 


lo.jD i 


o 
z 


1 Q OO 

io.zz 


Q 
O 


1 O 

iz.o ( 


o 
z 


ib.UU 


4 


TJT\ 0'7'7'7Q 

xlJJ Zl I lo 


1 Q /I Q 

lo.4o 





iz.oo 


c 



io.io 





1 c: OO 
10. oz 





XlJJ 0Z04:Z 


io.oo 


1 


<ii.OO 


/ 


1 Q /1 1 

io.41 


( 


iO.oU 


o 
o 


XlU / oooz 


lo.o ( 


y 


io.oo 


y 


io.oo 


iU 


1 /I o^ 


o 
o 


XlU 900/0 


1 O Q /I 

io.o4 


-1 -1 

ii 


io.4o 


-1 -1 

ii 






>io.UU 


-1 -1 

ii 


MU iUzUbo 


1 O To 


-1 -1 

ii 


1 o no 
io.Uo 


-1 -1 

ii 






1 Q OC 

io.oo 


-1 -1 

11 


XlU lUoyz/ 


















HD 110432 


13.19 


12 


13.25 


12 


12.20 


10 


14.30 


13 


HD 154368 


13.80 


2 


13.67'' 


3 


13.26 


2 


16.00'' 


8 


HD 167971 


13.53 


3 


13.73 


3 










HD 168076 


















HD 170740 


13.64 


9 


13.26 


9 


12.78 


9 






HD 185418 


13.12 


2 


13.10 


2 










HD 192639 


13.45 


2 


13.61 


2 










HD 197512 


















HD 199579 


13.36 


14 


13.01 


14 


12.09"^ 


15 






HD 203938 


13.61 


2 


13.68'' 


3 


13.31 


2 






HD 206267 


13.41^^ 


16 


13.02 


17 


12.91 


2 






HD 207198 


13.56 


2 


13.18 


17 


12.65 


2 


15.41'^ 


18 


HD 207538 


13.63 


2 






12.96 


2 


15.15'' 


19 


HD 210121 


13.48 


20 


12.78 


21 


13.09 


10 


15.48'' 


20 


HD 210839 


13.31 


22 


13.17 


22 


12.57 


2 


15.15"^ 


18 



''Emission measurement which should be considered only an approximation to the column density along 
the pencil-beam toward the disk of the star sampled by absorption measurements. 

''The referenced author(s) assumed h = 2 km 

'^The referenced author(s) assumed 6 = 1 km 

References. — (1) Snow et al 1994; (2) D. W. Welty, et al. 2002, in prep.; (3) Allen 1994; (4) 

Kaczmarczyk 2000; (5) Meyer & Roth 1991; (6) Joseph et al. 1986; (7) Cardelli et al 1990; (8) van 
Dishoeck et al 1991; (9) Gredel, van Dishoeck, & Black 1993; (10) Gredel, van Dishoeck, & Black 1991; 
(11) Gry et al. 1998; (12) Crawford 1995; (13) Paper II; (14) Jenniskens et al. 1992; (15) Joseph, Snow, 
& Scab 1989; (16) Fedcrman et al. 1994; (17) Chaffee & Dunham 1979; (18) Federman & Lambert 1988; 
(19) Dickman et al. 1983; (20) Welty & Fowler 1992; (21) de Vries & van Dishoeck 1988; (22) Crane, 
Lambert, & Sheffer 1995; 
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Table 3. Extinction parameters 



Star 


E{B - V) 


Ref 


Phot.'' 


Filters 


Ref 


Polar.'' 


Ref 


E.C.'= 


Ref 


Adopted 


Av 


BD +31° 643 


0.84 


1 


3.19 


HKL 


2 


3.75 


3 


3.46 


1 


3.19 


2.68 


HD 24534 


0.45 


4 








3.47 


5 


3.20 


4 


3.47 


1.56 


HD 27778 


0.38 


6 


2.65 


JK 


6 










2.65 


1.01 


HD 62542 


0.37 


7 


2.90 


JHKL 


7 


3.27 


8 


2.14 


9 


2.90 


1.07 


HD 73882 


0.72 


10 


3.16 


JHKL 


11 


3.51 


8 


2.93 


10 


3.16 


2.28 


HD 96675 


0.31 


12 


3.45 


JHKL 


12 


2.80 


13 


4.02 


12 


3.45 


1.07 


HD 102065 


0.17 


12 


4.24 


JHKL 


12 






3.36 


12 


4.24 


0.72 


HD 108927 


0.23 


12 


2.66 


JHK 


12 






3.73 


12 


2.66 


0.61 


HD 110432 


0.40 


14 


d 


JHKL 


11 


3.30 


15 






3.30 


1.32 


HD 154368 


0.82 


16 


3.02 


JHK 


17 






3.14 


16 


3.02 


2.48 


HD 167971 


1.04 


18 


3.30 


JHKL 


18 










3.30 


3.43 


HD 168076 


0.79 


19 


d 


JHKL 


19 


3.19 


20 


3.62 


21 


3.19 


2.86 


HD 170740 


0.48 


22 


2.61 


JHKL 


11 


3.08 


15 






2.61 


1.25 


HD 185418 


0.51 


23 












3.98 


10 


3.98 


2.03 


HD 192639 


0.66 


24 


2.83 


JHKLM 


25 










2.83 


1.87 


HD 197512 


0.33 


24 












2.56 


10 


2.56 


0.84 


HD 199579 


0.36 


24 


2.79 


KL 


26 






2.74 


10 


2.79 


1.00 


HD 203938 


0.72 


24 


3.04 


KLM 


26 






3.00 


10 


3.04 


2.19 


HD 206267 


0.52 


24 


2.63 


JHK 


27 










2.63 


1.37 


HD 207198 


0.62 


24 


2.20 


JK 


19 


2.30 


28 


2.66 


21 


2.20 


1.36 


HD 207538 


0.64 


24 








2.23 


8 






2.23 


1.43 


HD 210121 


0.38 


29 


2.10 


JHKL 


29 


2.13 


29 


2.01 


9 


2.10 


0.80 


HD 210839 


0.56 


24 


2.80 


JHKLM 


25 


2.86 


30 






2.80 


1.57 



''Derived from Eq. 1 using infrared photometry from the given reference 

''Derived from the wavelength of maximum polarization; Rv = 5.6Amax with Amax in /^m 

''Derived from the linear far-UV rise in the extinction curve 

''No reasonable solution could be obtained 

References. — (1) Snow et al. 1994; (2) Strom, Strom, & Carrasco 1974; (3) Andcrsson & Wannier 2000; (4) Snow 
et al. 1998; (5) Roche et al. 1997; (6) Kenyon, Dobrzycka, & Hartmann 1994; (7) Whittet et al. 1993; (8) Martin, 
Clayton, & Wolff 1999; (9) Wclty & Fowler 1992; (10) Fitzpatrick & Massa 1990; (11) Whittet & van Breda 1978; 
(12) Boulanger et al. 1994; (13) Whittet et al. 1994; (14) Paper II; (15) Serkowski, Mathewson, Ford 1975; (16) Snow 
et al. 1996; (17) The, Wesselius, & Jannsen 1986; (18) Leitherer et al. 1987; (19) Aiello et al. 1988; (20) Orsatti, 
Vega, & Marraco 2000; (21) Jenniskons & Greenberg 1993; (22) Clayton & Mathis 1988; (23) Sudzias & Bobinas 1984; 
(24) Garmany & Stencel 1992; (25) Castor & Simon 1983; (26) Sneden et al. 1978; (27) MorbideUi et al. 1997; (28) 
Anderson et al. 1996; (29) Larson, Whittet, & Hough 1996; (30) McDavid 2000 
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Tablc 4. Extinction curve parameters' 



Target 


(/xm~^) 


7 

(/Ltm~^) 


Cl 


C2 


C3 


C4 


Ref. 


BD +31° 643 


4.613 


1.28 


0.39 


0.51 


5.37 


0.58 


1 


HD 24534 




1.18 








0.552 


2 


HD 27778 
















HD 62542 


4.723 


1.174 


-2.416 


1.383 


1.949 


1.206 


3 


HD 73882 


4.576 


1.192 


-0.412 


0.788 


3.341 


0.540 


4 


HD 96675 


4.63 


1.10 


0.99 


0.35 


4.31 


0.58 


5 


HD 102065 


4.59 


1.19 


-1.50 


0.81 


8.00 


-0.21 


5 


HD 108927 


4.67 


1.05 


1.48 


0.44 


4.00 


0.65 


5 


HD 110432 
















HD 154368 


4.581 


1.15 


-0.01 


0.680 


4.419 


0.807 


6 


HD 167971 
















HD 168076 


4.595 


0.93 


0.57 


0.48 


2.85 


0.46 


1 


HD 170740 
















HD 185418 


4.579 


0.927 


1.266 


0.362 


3.941 


0.381 


4 


HD 192639 
















HD 197512 


4.585 


1.006 


-1.043 


1.021 


4.659 


0.438 


4 


HD 199579 


4.606 


0.997 


-0.725 


0.898 


2.923 


0.453 


4 


HD 203938 


4.589 


1.016 


0.087 


0.747 


3.647 


0.306 


4 


HD 206267 
















HD 207198 


4.596 


0.94 


-0.91 


0.95 


3.13 


0.66 


1 


HD 207538 
















HD 210121 


4.603 


1.085 


-2.493 


1.528 


2.405 


0.873 


3 


HD 210839 

















'^In the parameterization scheme of Fitzpatrick & Massa 1990 



References. — (1) Jenniskens &l Greenberg 1993; (2) Snow et al. 1998; (3) 
Welty & Fowler 1992; (4) Fitzpatrick & Massa 1990; (5) Boulanger et al. 1994; 
(6) Snow et al. 1996 
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Table 5. FUSE observations 



Target 


Date 




tint 

(ksec) 


S/N^ 


BD +81° 648 


1999 


Nov 13 


17 


37.5 


1.0 




2000 


Oct 26 


1 8 


84 7 


1.1 


HD 24584 


2000 


Sep 12 


8 


8.3 


24.3 


Hn 97778 


2000 


Oct 27 


4 


Q 7 


10 8 


HD 69549 


2000 


Jan 25 




1 1 


8 1 


HD 78882 


1999 


Oct 30 




25 5 


5 1 




2000 


Jan 24 


R 


1 1 q 


4 7 




2000 


Mar 19 


« 


18 6 


4 6 


HD Qfifi75 


2000 


Apr 10 


q 


10.2 


8 4 


HD 109065 


2000 


May 28 


a 


6 8 


7 8 


HD 108Q97 


2000 


May 27 


Li 


6 ^ 


1 1 9 


HD 1 1 0489 


2000 


Apr 4 




8 6 


98 ^ 


HD 154868 


2000 


Jul 14 


s 

O 


12.5 


2 8 


HD 1fi7Q71 


2000 


May 16 




Q 5 


2 


HD 168076 


2000 


May 16 


2 


6 6 


5 9 


HD 170740 


2001 


Apr 27 


r. 


9 q 


8 8 


HD 18^418 


2000 


Aug 10 




4 4 


1 4 q 


HD 192639 


2000 


Jun 12 


2 


4.8 


8.1 


HD 197512 


2001 


Jun 6 


3 


8.2 


10.9 


HD 199579 


2000 


Jul 19 


8 


4.3 


29.2 


HD 203938 


2000 


Jul 20 


4 


7.8 


3.5 


HD 206267 


2000 


Jul 21 


3 


4.9 


10.2 


HD 207198 


2000 


Jul 23 


3 


13.2 


11.2 


HD 207538 


1999 


Dec 8 


7 


7.7 


6.2 




2000 


Jul 21 


10 


11.2 


7.1 


HD 210121 


2000 


Jun 29 


5 


13.8 


2.9 


HD 210839 


2000 


Jul 22 


10 


6.1 


24.0 



^Number of integrations 
^Total integration time 

^Average per- pixel S/N for a 1 A region of the LiF 
lA spectrum near 1070 A, between the Lyman (3,0) and 
(2,0) bandhcads of H2. One resolution element corre- 
sponds to about 9 pixels. 



Table 6. Fitting test for HD 199579 



Segment Band Method Weights log A^(0) log iV(l) 

(cm (cm 



LiF 


lA 


(2,0) 


CURFIT 


None 


20.291 


20 


155 


LiF 


lA 


(2,0) 


AMOEBA 


None 


20.290 


20 


154 


LiF 


lA 


(2,0) 


CURFIT 


S/N 


20.286 


20 


154 


LiF 


lA 


(2,0) 


AMOEBA 


S/N 


20.285 


20 


153 


LiF 


lA 


(2,0) 


CURFIT 




20.308 


20 


157 


LiF 


lA 


(2,0) 


AMOEBA 




20.306 


20 


155 


LiF 


lA 


(4,0) 


CURFIT 


None 


20.319 


20 


106 


LiF 


lA 


(4,0) 


AMOEBA 


None 


20.319 


20 


105 


LiF 


lA 


(4,0) 


CURFIT 


S/N 


20.314 


20 


092 


LiF 


lA 


(4,0) 


AMOEBA 


S/N 


20.314 


20 


091 


LiF 


lA 


(4,0) 


CURFIT 


x' 


20.336 


20 


134 


LiF 


lA 


(4,0) 


AMOEBA 


x' 


20.335 


20 


132 


LiF 


2A 


(1,0) 


CURFIT 


None 


20.314 


20 


138 


LiF 


2A 


(1,0) 


AMOEBA 


None 


20.316 


20 


141 


LiF 


2A 


(1,0) 


CURFIT 


S/N 


20.323 


20 


144 


LiF 


2A 


(1,0) 


AMOEBA 


S/N 


20.321 


20 


141 


LiF 


2A 


(1,0) 


CURFIT 




20.297 


20 


123 


LiF 


2A 


(1,0) 


AMOEBA 


x' 


20.289 


20 


116 
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Table 7. Band-to-band variation 



Segment 


Band 










di^ 




LiF lA 


(2,0) 


-1 


44±0 


15 


-0 


39±0 


21 


LiF lA 


(4,0) 


+0 


52±0 


17 


-0 


40±0 


23 


LiF 2A 


(1,0) 


+0 


93±0 


16 


-0 


36±0 


25 


LiF 2B 


(4,0) 


-0 


29±0 


16 


+0 


40±0 


18 


SiC lA 


(2,0) 


-0 


12±0 


17 


+0 


24±0 


26 


SiC lA 


(4,0) 


-0 


30±0 


17 


+0 


38±0 


17 


SiC 2B 


(1,0) 


+0 


68±0 


22 


-0 


40±0 


20 


SiC 2B 


(2,0) 


-0 


14±0 


20 


+0 


55±0 


25 


SiC 2B 


(4,0) 


+0 


17±0 


17 


-0 


03±0 


24 


All 


(1,0) 


+0 


80±0 


13 


-0 


38±0 


16 


All 


(2,0) 


-0 


56±0 


13 


+0 


15±0 


14 


All 


(4,0) 


+0 


02±0 


09 


+0 


09±0 


11 


All 


All 





00±0 


08 





00±0 


08 



^Average normalized deviation from the mean 
and error of the mean {an-i / s/n) , for J=0 and 
J=l measurements 
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Table 8. Molecular and atomic hydrogen parameters 



Target 


Bands 


loe N(ll9) 


loe N(Q) 


loe N(l) 


T 

-'- kin 


loe iVfH 1) 


Ref 


iH2 








fern" 2 ) 






(K) 








BD +31° 643 


3 


21.09±0.19 


20.82±0.16 


20.76±0.24 


73±48 


21.38±0.30 


1 


0.51±0.26 


HD 


24534 


9 


20.92±0.04 


20.76±0.03 


20.42±0.06 


57± 4 


20.73±0.06 


2 


0.76±0.05 


HD 


27778 


9 


20.79±0.06 


20.64±0.05 


20.27±0.10 


55± 7 


20.98±0.30 


1 


0.56±0.20 


HD 


62542 


3 


20.81±0.21 


20.74±0.21 


19.98±0.14 


43±11 


20.93±0.30 


1 


0.60±0.28 


HD 


73882 


10 


21.11±0.08 


20.99±0.08 


20.50±0.07 


51± 6 


21.11±0.15 


3 


0.67±0.13 


HD 


96675 


9 


20.82±0.05 


20.63±0.04 


20.37±0.08 


61± 7 


20.66±0.30 


1 


0.74±0.18 


HD 


102065 


9 


20.50±0.06 


20.25±0.06 


20.15±0.06 


70± 9 


20.54±0.30 


1 


0.65±0.21 


HD 


108927 


9 


20.49±0.09 


20.30±0.09 


20.03±0.09 


60±10 


20.86±0.30 


1 


0.46±0.21 


HD 


110432 


9 


20.64±0.04 


20.40±0.03 


20.27±0.04 


68± 5 


20.85±0.15 


4 


0.55±0.11 


HD 


154368 


7 


21.16±0.07 


21.04±0.05 


20.54±0.15 


51± 8 


21.00±0.05 


5 


0.74±0.06 


HD 


167971 


4 


20.85±0.12 


20.64±0.10 


20.44±0.15 


64±17 


21.60±0.30 


6 


0.26±0.22 


HD 


168076 


9 


20.68±0.08 


20.44±0.08 


20.31±0.09 


68±13 


21.65±0.23 


2 


0.18±0.12 


HD 


170740 


7 


20.86±0.08 


20.60±0.05 


20.52±0.11 


70±13 


21.15±0.15 


2 


0.51±0.13 


HD 


185418 


9 


20.76±0.05 


20.34±0.04 


20.56±0.05 


101±14 


21.11±0.15 


3 


0.47±0.11 


HD 


192639 


9 


20.69±0.05 


20.28±0.05 


20.48±0.05 


98±15 


21.32±0.12 


2 


0.32±0.09 


HD 


197512 


9 


20.66±0.05 


20.27±0.05 


20.44±0.05 


94±14 


21.26±0.15 


3 


0.33±0.11 


HD 


199579 


9 


20.53±0.04 


20.28±0.03 


20.17±0.03 


70± 5 


21.04±0.11 


2 


0.38±0.09 


HD 


203938 


6 


21.00±0.06 


20.72±0.05 


20.68±0.08 


74± 9 


21.48±0.15 


3 


0.40±0.11 


HD 


206267 


9 


20.86±0.04 


20.64±0.03 


20.45±0.05 


65± 5 


21.30±0.15 


6 


0.42±0.11 


HD 


207198 


9 


20.83±0.04 


20.61±0.03 


20.44±0.04 


66± 5 


21.34±0.17 


2 


0.38±0.12 


HD 


207538 


18 


20.91±0.06 


20.64±0.07 


20.58±0.05 


73± 8 


21.34±0.12 


2 


0.43±0.10 


HD 


210121 


5 


20.75±0.12 


20.63±0.11 


20.13±0.15 


51±11 


20.63±0.15 


7 


0.73±0.11 


HD 


210839 


9 


20.84±0.04 


20.57±0.04 


20.50±0.04 


72± 6 


21.15±0.10 


2 


0.49±0.08 



References. — (1) Present work; iV(H I) = 5.8x10'^^ E{B - V) - 2N{R2); (2) Diplas & Savage 1994; (3) 
Fitzpatrick & Massa 1990; (4) Paper H; (5) Snow et al. 1996; (6) Present work; Lya profile fitting; (7) Welty & 
Fowler 1992; 21-cm emission measurement with possible systematic errors relative to the absorption measures 



-62- 



Table 9. Comparison with Copernicus results 



Target 


Satellite 


log iV(H2) 


log iV(0) 


log iV(l) 


Ref. 






(cin-2) 


(cm~^) 


(cin-2) 




HD 24534 


FUSE 


20.92±0.04 


20.76±0.03 


20.42±0.06 






Copernicus 


21.04±0.13 


20.78±0.11 


20.70±0.13 


1 


HD 199579 


FUSE 


20.53±0.04 


20.28±0.03 


20.17±0.17 






Copernicus 


20.36±0.18 


20.08±0.18 


20.04±0.18 


2 


HD 210839 


FUSE 


20.84±0.04 


20.57±0.04 


20.50±0.04 






Copernicus 


20.78 






3 



References. — (1) Mason et al. 1976; (2) Savage et al. 1977; (3) Bolilin 
et al. 1978 
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Table 10. CH component structure^ 



Target 


(km s~^) 


^ ^ couip 




/h2 


BD +31° 643 


3.6 


1 


1.00 


0.51 


HD 24534 


1.7 


2 


0.72 


0.76 


HD 27778 


1.3 


2 


0.53 


0.56 


HD 73882 


3.6 


1 


1.00 


0.67 


HD 110432 


0.3 


2 


0.68 


0.55 


HD 154368 


1.5 


4 


0.78 


0.74 


HD 167971 


1.5 


4 


0.47 


0.26 


HD 170740 


2.5 


1 


1.00 


0.51 


HD 185418 


1.7 


2 


0.59 


0.47 


HD 192639 


1.7 


2 


0.56 


0.32 


HD 199579 


1.7 


5 


0.64 


0.38 


HD 203938 


1.7 


3 


0.45 


0.40 


HD 206267 


1.3 


3 


0.50 


0.42 


HD 207198 


1.3 


3 


0.74 


0.38 


HD 207538 


2.0 


4 


0.39 


0.43 


HD 210121 


3.6 


1 


1.00 


0.73 


HD 210839 


0.6 


2 


0.62 


0.49 



^Data from D. E. Welty et al 2002, in preparation, 
except Crawford 1995 for HD 110432 and Crane et al. 
1995 for HD 210839 

Spectral resolution 
'^Number of components 

'^Fractional abundance of strongest component 
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Table 11. Separation into "translucent" and "diffuse" components^ 



Target 


E{B - V) 


Htot 


/h2 


Toi 


log JV(H I) 


log Ar(H2) 


log Af(0) 


log Ar(i) 






(cm-2) 




(K) 


(cm-2) 


(cm-2) 


(cm-2) 


(cm-2) 



BD+31 643 
HD 24534 
HD 27778 
HD 62542 
HD 73882 
HD 96675 
HD 102065 
HD 108927 
HD 110432 
HD 154368 
HD 167971 
HD 168076 
HD 170740 
HD 185418 
HD 192639 
HD 197512 
HD 199579 
HD 203938 
HD 206267 
HD 207198 
HD 207538 
HD 210121 
HD 210839 



0.20/0.64 
0.20/0.25 
0.15/0.23 
0.20/0.17 
0.35/0.37 
0.14/0.17 
0.05/0.12 
0.06/0.17 
0.08/0.32 
0.40/0.42 
0.14/0.90 
0.08/0.71 
0.12/0.36 
0.04/0.47 
0.04/0.62 
0.04/0.29 
0.06/0.30 
0.15/0.57 
0.14/0.38 
0.13/0.49 
0.13/0.51 
0.15/0.23 
0.11/0.45 



21.06/21.57 
21.06/21.02 
20.94/21.12 
21.06/20.99 
21.31/21.26 
20.91/20.99 
20.46/20.84 
20.54/21.00 
20.67/21.05 
21.37/21.20 
20.91/21.66 
20.67/21.70 
20.84/21.34 
20.37/21.34 
20.37/21.45 
20.37/21.40 
20.54/21.15 
20.94/21.62 
20.91/21.42 
20.88/21.44 
20.88/21.49 
20.94/20.83 
20.80/21.33 



0.90/0.38 
0.90/0.60 
0.90/0.35 
0.90/0.25 
0.90/0.41 
0.90/0.61 
0.90/0.54 
0.90/0.30 
0.90/0.41 
0.90/0.51 
0.90/0.15 
0.90/0.11 
0.90/0.38 
0.90/0.43 
0.90/0.27 
0.90/0.28 
0.90/0.26 
0.90/0.29 
0.90/0.27 
0.90/0.25 
0.90/0.31 
0.90/0.50 
0.90/0.37 



30/187 
30/185 
30/158 
30/107 
30/163 
30/160 
30/140 
30/119 
30/191 
30/196 
30/166 
30/142 
30/175 
30/158 
30/172 
30/165 
30/192 
30/165 
30/170 
30/174 
30/177 
30/138 
30/170 



20.06/21.36 
20.06/20.62 
19.94/20.94 
20.06/20.87 
20.31/21.04 
19.91/20.58 
19.46/20.50 
19.54/20.84 
19.67/20.82 
20.37/20.89 
19.91/21.59 
19.67/21.65 
19.84/21.13 
19.37/21.10 
19.37/21.32 
19.37/21.25 
19.54/21.03 
19.94/21.47 
19.91/21.28 
19.88/21.32 
19.88/21.32 
19.94/20.53 
19.80/21.13 



20.72/20.85 
20.72/20.50 
20.59/20.36 
20.72/20.09 
20.96/20.58 
20.56/20.47 
20.12/20.28 
20.19/20.18 
20.32/20.36 
21.02/20.60 
20.56/20.54 
20.32/20.43 
20.50/20.62 
20.02/20.68 
20.02/20.59 
20.02/20.55 
20.19/20.26 
20.59/20.79 
20.56/20.55 
20.53/20.45 
20.53/20.68 
20.59/20.23 
20.46/20.60 



20.70/20.19 
20.70/19.84 
20.58/19.75 
20.70/19.63 
20.95/19.96 
20.55/19.86 
20.10/19.74 
20.18/19.68 
20.31/19.69 
21.01/19.92 
20.55/19.91 
20.31/19.86 
20.48/19.97 
20.01/20.07 
20.01/19.95 
20.01/19.93 
20.18/19.59 
20.58/20.16 
20.55/19.91 
20.52/19.89 
20.52/20.03 
20.58/19.67 
20.44/19.97 



19.19/20.75 
19.19/20.39 
19.07/20.24 
19.19/19.90 
19.43/20.46 
19.04/20.35 
18.59/20.14 
18.67/20.01 
18.79/20.26 
19.49/20.50 
19.04/20.42 
18.79/20.30 
18.97/20.51 
18.49/20.56 
18.49/20.48 
18.49/20.44 
18.67/20.16 
19.07/20.67 
19.04/20.43 
19.00/20.42 
19.00/20.57 
19.07/20.09 
18.93/20.49 



''In each pair of values, the value for the "translucent" component is given first. 



